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prepared for direct photographic reproduction except for the lettering which should be 
inserted in pencil. Legends should be given in the manuscript indicating where in the 
text the figure should appear. Blocks are retained by the Society for 10 years; unless the 
author requires them before the end of this period they are then destroyed. 


5. Tables—These should be arranged so that they can be printed upright on the 
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6. Proofs.—Costs of alteration exceeding 5 per cent of composition must be borne 
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simplification of symbols and equations consistent with clarity. 
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MONTHLY NOTICES 
OF THE 


ROYAL ASTRONOMICAL SOCIETY 
Vol. 110 No. 3 


MEETING OF 1950 MARCH 10 
Professor W. M. Smart, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Reginald Gordon Andrews, The Old Brick House, Town Littleworth, 

_ Cooksbridge, Lewes, Sussex (proposed by F. M. Holborn); 

Henry Reginald Calvert, M.A., B.Sc., D.Phil., The Science Museum, 
South Kensington, S.W.7 (proposed by R. d’E. Atkinson) ; 

Herbert Chatley, D.Sc., M.I.C.E., A.Inst.P., 4 Belgrave Road, Grosvenor, 
Bath (proposed by A. T. Doodson) ; 

Arthur D. Code, Leander McCormick Observatory, University of Virginia, 
Charlottesville, Virginia, U.S.A. (proposed by S. Chandrasekhar) ; 

Sidney Charles Bartholomew Gascoigne, Commonwealth Observatory, 
Mount Stromlo, Canberra, Australia (proposed by R. v. d. R. Woolley); 

John Herbert Hewitt, 60 Farquhar Road, Edgbaston, Birmingham (proposed 
by M. C. Johnson); 

Arthur Robert Hogg, Commonwealth Observatory, Mount Stromlo, 
Canberra, Australia (proposed by R. v. d. R. Woolley) ; 

John Alexander Jenkins, Glaselune, Balerno, Midlothian, Scotland (proposed 
by W. M. H. Greaves) ; 

Ralph Leslie Maughan, M.Sc., F.Inst.P., Constantine Technical College, 
Middlesbrough, Yorkshire (proposed by F. R. Curtis) ; 

Alexander Charles Robbins, 64 Eastbury Road, Northwood, Middlesex 
(proposed by R. C. Wood); and 

Leonard Stanley Theodore Symms, 21 Shrewsbury Lane, Shooter’s Hill, 
London, S.E.18 (proposed by R. d’E. Atkinson). 


One hundred and twelve presents were announced as having been received 
since the last meeting, including :— 


Harlow Shapley, Galactic and Extra-Galactic Studies. First Collection 
(presented by the author) ; 
Peter Doig, A Concise History of Astronomy (presented by the author) ; 
Harold Jeffreys, Earthquakes and Mountains (presented by the author) ; 
érard de Vaucouleurs, The Planet Mars (presented by the author); and 
P. A. Schilpp, Ed., Albert Einstein: Philosopher Scientist. Volume VII of 
the Library of Living Philosophers (presented anonymously). 
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MEETING OF 1950 APRIL 14 


Professor W. M. Smart, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Alan Hugh Cook, Department of Geodesy and Geophysics, Downing Place, 
Cambridge (proposed by B. C. Browne) ; 

Sidney Guy Hacker, Department of Mathematics, State College of Washing- 
ton, Pullman, Washington, U.S.A. (proposed by H. N. Russell) ; 

John Meadows Jackson, B.Sc., Ph.D., University College, Dundee (proposed 
by Sir Edmund Whittaker) ; 

Gerry Marlowe, The Stables, Burford Drive, Whalley Range, Manchester, 16 
(proposed by J. C. Farrer); and 

Stanley Keith Runcorn, Gonville and Caius College, Cambridge (proposed 
by R. I. B. Cooper). 


One hundred and fourteen presents were announced as having been received 
since the last meeting, including :-— 


Angus Armitage, A Century of Astronomy (presented by the author) ; 

Oliver Justin Lee, Measuring our Universe (presented by the author); and 

An original photograph of Comet 19481 taken during the Eclipse (presented 
by R. O. Hennings). 


The President announced that the Council had gratefully accepted a legacy, 
amounting to approximately £750, from the late Mrs M. A. Nadarov, bequeathed 
to the Society in the following terms :— 

“To the Royal Astronomical Society, Burlington House, London, for the 


encouragement and promotion of the study of Sunspots, and to be known as the 
Victor Nadarov Fund.” 
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A STUDY OF THE ERUPTIVE PROMINENCE OF 
1948 SEPTEMBER 27 


Helen W. Dodson and Edwin B. Weston 


(Communicated by Robert R. McMath) 


(Received 1950 March 10) 
Summary 

Motions in an eruptive prominence are discussed. The radial velocities 

of the prominence gases are found to be consistent with the cross-motions 

shown by the spectroheliograms. This agreement is strong evidence in 


favour of the hypothesis of moving material. Motions of descending nodules 
and of an unusual ascending spike are included in the study. 


1. Introduction 


On 1948 September 26 and 27, a large prominence was observed to arch 
across the solar equator on the east limb of the Sun. A unipolar spot (Mount 
Wilson No. 9438), with field strength 1100 gauss, was located near its northern 
extremity. There was no spot near the southern extremity. Throughout 
September 26, faint knots or streamers passed from the southern to the northern 
branches of the prominence. This type of activity continued until 15" on 
September 27 when the entire bright southern branch of the prominence began 
to move upward and along the arch. As this bright mass advanced it disturbed 
the stability of the arch and the entire phenomenon became an eruptive 
prominence. Hydrogen spectroheliograms secured with the Tower Telescope 
of the McMath-Hulbert Observatory and spectra recorded by the Stone 
Spectroheliograph permit a detailed study of the prominence and provide evidence 
in favour of the hypothesis of moving material within the prominence. A well- : 
defined ascending spike was a secondary feature of considerable interest. 

Intercomparison of radial velocities and motions in the plane of the sky, 
coupled with the data secured on the two successive days, indicate that on 
September 27 the northern tip of the prominence extended over the disk while 
the southern extremity emanated from a point on or behind the limb. Motion 
from south to north along the arch of the prominence would thus contain a 
component directed towards the observer. 


2. Summary of Principal Activity 

(a) Between 15" 40™ and 17" 30™ the bright mass moved northward, upward 
and towards the observer with slowly increasing velocity. 

(b) During this period, at least part of the bright mass was in rotation about 
an axis that was more nearly radial than tangential to the solar surface. As 
viewed from the centre of the Sun the direction of rotation was counter- 
clockwise. The rotation is apparent when the film is viewed as a cinema and it 
is confirmed by the positive and negative radial velocities in the spectra. 

(c) As the bright mass advanced, it disturbed the relatively stable arch or 
streamer that had been observed since September 26. It should be noted that 
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prior to 17" 00™ there was no significant change in the position of the northern 
half of the arch even though the southern half had been completely modified. 

(d) Between 16°23™ and 17508" there developed within the prominence a 
conspicuous pattern which was concave upward. As the bright mass moved, the 
pattern moved ; its orientation was fixed with respect to the moving prominence, 
not with respect to the chromosphere or photosphere. Fig. 1 (Plate 3) shows 
this pattern at 16" 29™ and 16556™, 

(e) Material continually dropped back to the chromosphere, principally along 
the northern and southern branches of the original arch, but in part along new 
paths emanating from the bright mass. The descending material along the 
northern branch was more intense than that along the southern branch. This is. 
clearly shown in the series of spectra for 17" 28™, 

(f) In the final stages of the activity (17"25™ to 19"00™) the “bright mass” 
and the remnant of the stable arch completely merged. ‘The configuration 
expanded, and as it rose in height it continued to feed material back to the disk. 
The prominence was measured to a height of 500,000 km. above the solar surface. 
After 18" 35™ the spectroheliograms showed only descending material. 


3. Velocities and Accelerations in Principal Activity 


Although the phenomenon was complex and many of the features ephemera: 
or diffuse, it was possible to measure with a certain confidence the x and y 
coordinates of a limited number of regions. In the system of rectangular 
coordinates in which the measurements were made, the origin was on the solar 
limb in latitude 3° south, the x-axis was tangential to the solar disk at this point, 
and the y-axis passed through the origin in the direction of the radius of the 
apparent solar disk.* ‘The measures of x and y are plotted versus time in Fig. 2. 
Velocities in the x and y directions have been determined from the slopes of 
curves drawn through the mean of the plotted points. It should be noted that in 
this study the measurement of x was considerably less certain than that of y. 
Velocities in the z, or radial, direction have been determined from spectro- 
heliograms taken with the Stone Spectroheliograph.t The radial velocities 
are based on displacements of Ha emission of the prominence measured with 
respect to the dark Ha line of the scattered light from the adjacent solar limb. 
The effect of solar rotation on the radial velocities is thus eliminated. The 
spectra in Fig. 1 (Plate 3) are oriented in such a manner that negative radial 
velocities are represented by displacements downward; positive radial velocities 
result in displacements upward. Radial velocities are not available for the highest 
parts of the prominence nor for the very faint portions within it. In the former 
case the regions were beyond the field of the Stone instrument; in the latter case 
the spectra were too weak for satisfactory registration. 

During the period 15" 40™ to 17" 30™ the velocities in the x, y and 2 directions 
were of the same order of magnitude (10 to 30 km./sec.). The values are recorded 
in Table I. The magnitude and direction of the velocities in the xy-plane (the 
plane of the sky) are shown graphically in Fig. 3(a) against the background of the 
prominence. The tracks of the arrows in this diagram trace the trajectories of 

* A more detailed description of this coordinate system is given in Publ. Obs. Univ. Mich., 8, 
123, 1941. 


t Detailed descriptions of the instrument and the method of measurement are given in Publ. 
Obs. Univ. Mich., 8, 57, 1940. 
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Fic. 1.—H« spectroheliograms and spectra for the prominence of 1948 September 27. 
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the measured features. Concurrent radial velocities in km./sec. are indicated 
by the adjacent numbers. 

Between 17" 30™ and 18" 30, velocities in the y, or height, direction increased 
rapidly. For all ascending features the last measurable spectroheliograms 
indicate outward velocities between 150 and 200km./sec. The similarity of the 
y versus time curves for the various features should be noted. 


Fic. 2.—Plot of x and y for all measured features in the prominence. 


The accelerations in the xy-plane represented by these velocities are likewise 
recorded in Table I and shown graphically in Fig. 3(4). It should be emphasized 
that accelerations are measures of second differences and of necessity carry a 
high degree of uncertainty. There are, nevertheless, certain aspects that seem 
clearly indicated. Throughout the entire activity, accelerations in the x direction 
were small; the same is indicated by all available data for the z-plane. 
Accelerations in the y, or height, direction dominated the activity and these began 
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Positions, Velocities and Accelerations for the Measured Features of the Prominence of 1948 September 27 


Feature 
measured Sept. 27 


Position 


(km. 10%) 


Velocity 
(km./sec.) 


Acceleration 
(km./sec.? x 10—*) 


h m 


Top point of 15 45 
ascending 50 
bright mass 55 
16 00 

A 05 

10 

15 


17 00 


x 


Vez Vy V; 
Principal Activity 


14 12 —10 
14 —II 
15 12 —I! 
15 13. —12 
16 13. —%3 
16 13° —13 
16 13. —14 
16 13. —I5 
17 14 —I5 
17 14 —16 
17 14 —17 
17 14 —18 
17 14 —19 
17 16 —20 
17 18 —2I 
17 19 
17 20 —22 
17 21 —23 
18 22 —23 
18 23 —24 
19 28 —25 


20 31 —26 


—28 
14 42 —30 
—40 


—28 
23 4§ 
—40 


—28 
14 34. —30 


—50 


No FWNDN 


A\/A3+ A} 


—2 
—3 
—3 
—3 
=—3 


Vol. 110 


—7o 156 18 21 06 08 

—66 160 18 22 of o8 

—62 164 19 22 06 

—57 168 20 23 06 

—52 172 21 24 06 

Ber —48 176 21 24 06 o8 

—44 180 21 25 of 

20 —38 184 21 26 of o8 

25 —33 188 22 27 06 08 

30 —28 192 22 27 06 o8 

35 —23 196 28 of 

40° —17 200 22 29 06 I 

45 —12 204 22 29 

50 (—8) 208 23 31 

55 (—2) 213 25 32 

(4) 219 25 33 O7 

05 (8) 225 26 34 

10 (13) 231 27 35 

15 19 237 28 37 «08 

20 24 244 29 38 «008 

25 30 252 34428 : 

Ascending 

| nodules 17 40 62 276 44 52 

| 50 70 302 44 60 

45 54 234 5° 59 

50° 62 248 50064 

D 17 40 80 254 37.46 

45 84 264 37.48 

5° 88 274 14 34 —-40 37 54 

E 18 00 T7° 292 to 722 to 68 

F 18 00 116 252 22 34. —50 «664 29 29 

02°5 120 258 22 40 

128 270 8 22 530 55 55 

10 132 278 22 60... 60 

12°5 135 288 22 FB ose 82 

15 138 300 22 103... 103 

20 146 336 22 188... 188 
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Tasie I—Continued 
Positions, Velocities and Accelerations for the Measured Features of the Prominence of 1948 September 27 


Feature pi Position Velocity Acceleration 
measured pa (km. x 10°) (km./sec.) (km./sec.? x 107%) 


h m X Ve :¥Vy Ve Valen fy 
Ascending Principal Activity 
nodule 
G 18 00 308 6 68 (65) 


115 180 
Descending 19 —121 86 
material 19 —I21 t 
after 19 —I2I 
eruption 5 19 —I121 
H 5 19 —121 
19 —I2I 


Secondary Activity 

Descending 
ejected from —43 (—25) 48 (5 —66 
ascending 27° —sr (—25) — 66 
bright mass 30 —59 (—25) 60 ; — 66 
R 32° 5 —66 
35 - 5 —86 (—25) 86 —66 
—94 (—25) 95 —66 
40° - —104 (—25) 104 —66 
42°5 —112 (—25) 115 —66 


Ascending 17 10 I4 116 —46 117 189 —30 
spike 12°5 — 21 154 —51 156 189 —30 
Ss 15 30 198 —56 201 189 —30 

iy 34 201 —61 204 189 —30 

20 50 198 204 214 189 —30 


( ) Uncertain value. 

* Measures cover too short a time interval for determination of acceleration. 

+ No significant change in the value of X, but measures are too uncertain to permit meaningful 
determinations of velocity and acceleration. 

t The radial velocity for this object is positive but its change in value with time is unknown. 


to increase in magnitude at about 17"00™. Within the range of the measurements 
the accelerations in the y direction either remained constant or increased with 
increasing height. 

After 18" 30™, the principal feature was a relatively bright descending nodule 
(Feature H). This appeared to move in the xy-plane with the large constant 
velocity of 123km./sec. Its radial velocity at 188 30™ was + 86km./sec. 


4. Secondary Activity 
(1) Descending nodules ejected from ascending bright mass.—At about 


16 15™ a series of knots (Feature R) was ejected from one of the rotating “‘ veins ” 
of the ascending mass. As the knots descended they formed a streamer and 


123 150 ° ° 
: ° ° 
° ° 
23 ° ° 
73 i 
100 
72 
198 
198 
198 
198 
198 
pie 
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entered the chromosphere near latitude 10° south. The individual knots followed 
curved trajectories in the xy-plane, but their locus at any one instant was essentially 
linear. The feature is clearly shown in the spectroheliogram for 16% 29™ in 
Fig. 1. At 16%15™ the point from which the nodules emanated showed the 
presence of gases with radial velocities of the order of —8okm./sec. It is not 
possible to determine the variation of V, with time for the individual nodules, 
but the data are sufficient to show that the radial velocities remained negative 
and between 16" 25™ and 16" 45™ were of the order of — 25 km. sec. 

The x and y measures for the first knot ejected are shown as Feature (R) in 
Table I and Figs. 2 and 3. The velocities increased with time and the space 
accelerations showed a relatively larger component at right angles to the 
trajectory than was the case for the ascending material. 


T 


. 


Fic. 3.—(a) Velocities: magnitudes and direction of velocities in the xy-plane are shown by 
arrows; velocities in the z or radial direction are indicated by numerals. (b) Accelerations : 
magnitude and direction of accelerations in the xy-plane are shown by arrows. 


(2) Ascending spike.—Another secondary ejection (Feature S) took place 
at 17"10™. Material was apparently shot some 100,000 km. upward from a 
relatively immobile portion of the bright mass. Its space velocity was very great 
and increased from 125 to 214km./sec. The acceleration was very large, almost 
o-2km./sec.”, and was directed outward. Its form was that of an ascending 
bright spike and in this respect it could be said to resemble the surge-type 
prominence. However, the properties of this spike and those of the surge-type 
prominence are so different that they cannot be identified with the same 
activating mechanism :— 


(a) The spike had its base in an elevated prominence, not in the chromosphere 
or photosphere. 

(6) There was no evidence of the return of the material along the path of 
ejection, with consequent diminution in apparent height. 

(c) The ascending motion showed an acceleration upward rather than 
downward, as is the case for the surge-type prominence. 
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5. Evidence for Moving Material 


Intercomparisons of the geometry of the prominence, the velocities in the 
x, y and z directions and their changes with time provide the following data 
which are consistent with the hypothesis of moving material : (1) Between 
15" 45™ and 17" 30™ the velocities in the x, y and z directions were of the same 
-order of magnitude and all increased with time. (2) As the bright mass moved 
upward over the southern branch of the arch, its radial velocity was negative. 
After it had risen and material had begun to descend to the south along 
apparently the same trajectory, the radial velocity was positive. This change in 
the sign of the radial velocity can be seen clearly by comparison of the spectra 
for 16" 02™ and 17" 28™ in Fig. 1 (Plate 3). In the former the bright mass was 
ascending the southern branch of the prominence and at a height of 136,000 km. 
(well within the body of the feature, not its upper tip) the radial velocity was of 
the order of —30km./sec. Inthe later scene, at a corresponding position on the 
trajectory, the radial velocity was of the order of +70km.)sec. These later spectra 
also show material descending on the northern branch of the prominence with a 
negative radial velocity of about 50km./sec. ‘These changes in the sign of the 
radial velocities are exactly those that would be expected if actual material were 
moving in the manner here described along an arching prominence whose 
southern tip was behind the limb of the Sun and whose northern extremity was 
already projected on the solar disk. 
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MONOCHROMATIC MAGNITUDES OBSERVED WITH A 
COARSE GRATING 


R. v.d. R. Woolley and K. Gottlieb 
(Received 1950 January 16) 


Summary 


Very low dispersion spectra obtained with a coarse grating in front of a 
30-inch reflector were used to determine magnitudes at 5980 A. and 4480 A. 
From internal evidence, and from comparisons (rather few in number) with 
Hogg and Mason’s photoelectric work, the standard error of the mean 
magnitudes derived is about +o™-o5. Comparison with Dr Gascoigne’s 
work indicates that the standard error of the gradients is +0-°08. The 
method is capable of extension to the tenth magnitude without unduly long 
exposure time. 


1. The Reynolds 30-inch reflector of the Commonwealth Observatory at 
Mount Stromlo is equipped with a coarse grating. As the grating has 1o wires 
in 3cm., and as the focal length of the Newtonian arrangement of the telescope 
is 120 inches, the grating produces little spectra at this focus with a dispersion 
of g800 A.mm. This dispersion is of course very tiny, but the little spectra 
have certain advantages over objective prism spectra made with a refracting 
telescope, since the scale is not cramped in the red relative to the blue and there 
is no “secondary spectrum’’ due to chromatic aberration in the telescope. 
: Also, the wave-length of any part of the spectrum is fixed by its distance from 
; the central image. The resolving power in the spectra is limited by plate grain 
and by the performance of the telescope, and depends on the state of the figure of 
; the mirror, the seeing and the guiding. Under good conditions, the resolution 
is such that the image of a truly monochromatic source would fall off to half the 
central intensity at about 250A. distant on either side of the centre. Notwith- 
standing the tiny scale and feeble resolution it was decided to measure some 
grating spectra with slits sampling about 250A. of the photographic images, in 
order to see whether magnitudes or gradients could be determined with a 
reasonably small error; because if the method were at all successful it would 
give a means of measuring the quantities concerned in stars down to the ninth 
or tenth magnitude without unduly prolonged exposures. 

It was found at the outset that the consistency of the results was improved by | 
trailing the spectra in the direction perpendicular to dispersion. The grating 
was therefore set to disperse perpendicular to the direction of daily motion and 
the telescope was trailed by driving it at three-quarters of the diurnal rate. Since 
the speed at which the image moved on the plate depended on the declination of 
the star it was necessary to perform experiments with comparison spectra trailing 
at ditferent speeds to determine the reciprocity failure of the photographic plates. 

The spectra of various stars were brought to comparable densities by spinning 
appropriate sectors just inside the focus, a set of sectors cutting down by 
5 magnitudes, 44 magnitudes and so on in } magnitude steps down to 4 magnitude 
being used. At least two sectors (counting no sector at all as the “0 sector”). 
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were spun in front of each star on each plate. Each photographic plate used. 
was calibrated in an auxiliary spectrograph, but the method depends for its. 
accuracy on the sectors used in the stellar exposures. 

The spectra appear back to back symmetrically disposed about a central. 
trail. They were measured in a projection photometer carrying a double slit in 
front of the photocell, the distance between the slits being twice the wave-length 
required on the projected scale of the spectra. The use of a double slit makes 
setting easier and more reliable than is the case with a single slit. All the plates. 
were measured with slit pairs set at 5980 a. and 4480 A., and about a quarter of the 
plates were also measured with a pair of slits set at 5780 a. Since the measures. 
with the intermediate slit pair interpolated well between the two extremes, it was 
decided not to measure all the plates in the extra wave-length. The widths of all 
the slits used were about 250A. on the projected scale of the spectra. 

No attempt was made to determine atmospheric absorption. Instead, stars 
were observed with 1-4>secz>1-2 as far as possible, and corrections of 
o™-3 secz in the blue and o™-1 sec z in the red applied. 

The results were tested by comparing the magnitudes of eight stars with 
photoelectric observations (unpublished) by Mr Hogg and Miss Mason at 4550 A. 
The blue coarse grating magnitudes uncorrected for the nominal wave-length 
difference of 70 a. show a standard error of + 0™-06 tested against these; when a 
correction for the wave-length difference is applied, the standard error falls to 
+0™-046. Since the spectral types range from Bo to K2, the nominal wave- 
length of the grating magnitudes cannot be far in error. In addition to this test, 
the observed gradients of eight stars were compared with (unpublished) 
determinations by Dr Gascoigne, the standard deviation in the gradient ¢ being 
+0°08. (Six stars were common to the two lists.) 

The exposure time depended only on the length of trail required, i.e. with a 
fixed trailing speed of the polar axis, on the declination of the star. Stars of 
— 60° declination were trailed for 20 seconds. Sinve the method at present reaches 
to the sixth magnitude, a modification which repeats the trails will go well down 
to the ninth magnitude with 10-minute exposures.* 

The least satisfactory features of the method are reciprocity failure and 
waste of light. The dependence of trailing speed upon declination necessitates 
a separate experiment to determine reciprocity failure, and it, is proposed to 
eliminate this feature by driving the plateholder backwards and forwards relative 
to the telescope, which will be driven at normal speed. This modification will 
admit longer effective exposure times, and extend the range of magnitudes 
downwards from the sixth. (The apparatus will be a one-dimensional 
“‘ Schraffierkasette ”, and only one dimension is available if the other is used 
for dispersion.) 

Light is wasted by the grating; only half the starlight ever gets past the wires, 
and of this two-thirds goes into the central image, so that less than a tenth gets 
into each of the two spectra. It is intended to remove the grating and substitute 
a narrow-angle reflecting prism at the Newtonian flat (placed between a pair of 
singlet lenses of opposite power to make the light parallel and focus it again). 
This modification will either give more dispersion or cut down exposures on faint 
stars; but it will sacrifice the linear dispersion. 


* Note added 1950 April 14.—The star H.D. 25224, magnitudes 9°27 red, 10°14 blue, has been 
measured using 15-minute exposures. 
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2. Each sector used in cutting down starlight is mounted permanently on the 
spindle of a small 6-volt D.C. motor which is carried on a substantial aluminium 
slide which may be pushed into a keyway in front of the plateholder. Electric 
contact is established only when the slider is in the correct position. Two types 
of plateholder were used in the series, taking quarter plates and quarter plates 
cut in half parallel to the longer side respectively. In the later plates in the series 
quarter plates cut in half were used ; on these twelve stellar spectra were exposed, 
at least one of which was a repeat of an earlier exposure. In addition four 
comparison spectra were exposed. Such a plate takes about an hour to secure, 
more than half of the time being spent on setting the telescope. 

The plates were Panatomic X developed for 4 minutes at 65 deg. F. in Kodak 
D26 developer. 

The calibrating spectrograph is illustrated in Fig. 1. Four pinholes were 
mounted in a slide A at the focus of a collimator. ‘This slide was moved 
horizontally by a worm and nut driven by a 50-cycle synchronous motor and 
suitable reduction gear. Parallel light from the collimator lens C passed onto a 
photograph D of the coarse grating. After passing through the photographic 
grating the light fell onto a small camera lens E having the same focal ratio as the 
Reynolds telescope, so that the spectra F of the pinholes were similar to those of 
the actual stars. 

In front of the pinholes A a step sector B is spun by a 6-volt motor similar 
to those used for the star sectors. This step sector runs in half magnitude steps 
from }™ to 3}™ and is mounted in a slide so that any four steps desired can be 
placed in front of the pinholes. The whole is illuminated by an automobile tail- 
light behind a blue filter and diffusing screen. In order to eliminate linear’errors 
in the illumination, comparison spectra were secured in pairs such as :— 

Stepp 41 1 2 
It was found best to adjust the distance of the lamp from the diffusing screen so 
that a hole with no sector in front of it gave a slightly overexposed trail, and then 
to make two pairs of trails with o™, #™, 1™, 13™ and 1™, 1}™, 2™, 2}™ sectors in 
front of the holes. There are then six points on the calibration curve for the ~ 
plate and the usable range of densities is well covered. 

The heights of the steps were equal to the distances between consecutive 
holes, and greater than the distance trailed by a comfortable margin. 

Intermittency.—Nearly all the spectra were secured with intermittent 
exposures, the only exceptions being the 0™ spectra in comparison images and 
some faint stars on which the sectors used were 0™ and 4” (or 0 and 1). No 
systematic difference was found.* 

Reciprocity failure.—Four worms and nuts, with pitches of J, inch, inch, 
j;in. and finch were employed to drive the pinhole slide A. These trailed the 
comparison spectra F at the same speed on the plate in mm./sec. as the telescope 
trailed stars of declination 82° 31’, 74° 48’, 58°22’ and 29°01’. Considering 
any particular point on the trail, the time integral of the amount of light received 


* The running of the motor depends on the state of the bearings and also on the position (relative 
to gravity) in which it is held. When running well the motor attains 3500 r.p.m., and when the 
rate falls so low as 2500 r.p.m., the motor labours audibly. When this occurred the motor was 
withdrawn from service next day and repaired or replaced. All sectors except one (5 mag.) carried 
two cuts so that the intermittency was, in these cases, doubled 
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is inversely proportional to the speed of the image, and the brightness of the 
images integrated in this way are 


+0m-76, +1™-51, + 2™-06, 


referred to the slowest moving (+ inch pitch) as zero. 

Seven plates were taken, on each of which were 64 spectra, 16 with each pitch. 
Each group of 16 spectra had four arrangements of the sectors and pinholes, 
and could be used in itself to calibrate the plate. Selecting any particular pitch 
to calibrate the plate, each set of 16 spectra taken with any other pitch gives 
(16 values of) the apparent brightness of the sources. All plates were measured 
in the red (6000 A.) and blue (45004.) and the mean results were as follows : 


Reciprocity Failure 
Pitch of Screw finch } inch 
Brightness with zero recipro- 
city failure o™00 76 2-06 
Apparent brightness, red 0-00 0°54 I-12 1-57 
Apparent brightness, blue 0-00 0°54 I-Il 1-59 


If m, is the brightness with zero reciprocity failure in magnitudes and m, 
ithe apparent brightness, the mean of the red and blue taken together is well 
represented by 


My = 0°75 


If the plates are taken individually and the best straight line adopted for each 
plate, the result is 


My==(0°77 + 0-03)m,, 


which indicates the reliability of the result. We adopt m,=0-76m, and then the 
correction to the apparent magnitude of a star of declination 4 is 


Am = 0-76 x 2:5 log sec3. 


Closing error.—The average magnitude difference between two exposures 
-of the same star taken on the same plate (usually about half an hour apart) is 
o™-07, from 75 such comparisons. This includes of course possible variations 
in the sky during the time taken to secure a plate, but the contributions from 
atmospheric reddening as distinct from thin cloud are not large, as the average 
failure to close is 0™-08, in the blue and 0™-06, in the red. If the errors were 
entirely atmospheric there would be more colour difference. 

3. The stars « TrAu, 8 TrAu, yTrAu, B Pavo, yPavo and 5 Pavo were 
taken as reference stars and intercompared frequently; the following relative 
magnitudes were adopted :— 


Blue Red 


TrAu 0-00 
BTrAu — 0:05 
y TrAu — 0°23 
B Pavo +0°39 
y Pavo +151 
5 Pavo +1-01 
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-and all the remaining magnitudes have been determined with reference to one 
of these stars. The results for 34 stars obtained from 44 plates exposed between 
1949 July 7 and 1949 September 17 are shown in Table I. The qualification 
for a star appearing in this table is its measurement on at least three platés. The 
zeros of the magnitude scales have been adjusted so as to make the average 
difference between the magnitudes of the four Ao stars observed and the 
magnitudes of these stars in the Harvard Mimeograph catalogue zero in both 
colours. No standard error is shown for «TrAu, and the standard errors 
shown for the remaining five reference stars are the standard errors of their 
magnitudes relative to « Tr Au. 

The table shows the results of 182 determinations. The sums of all the 
squares of the residuals of the individual determinations of a particular star 
from the adopted value for the star are 20,988 in the blue and 26,042 in the red 
{the unit being a hundredth of a magnitude squared), so that the average value 
-of Xv?/n is found by dividing those numbers by 182. The results are 115 (blue) 
and 143 (red). The standard errors of the means of observations are therefore 


and V143/(n—1) 


respectively, the unit being a hundredth of a magnitude. For example, with 
n=6, the average standard errors are 


+o0™-048 and +0"-053 


in the blue and red respectively. 

Monochromatic magnitudes for stars in the Hyades and in the region of 
Praesepe have been determined by Ramberg (Stockholms Observatoriums Annaler, 
13, No. 9, 1941). He used an objective prism attached to the 40-cm. refractor 
at Stockholm. On account of chromatic aberration in the object glass, only a 
portion of the spectrum was in good focus, and the magnitudes were determined 
at 42154. for spectral classes B to F and at 4360a. for classes G to M. The 
equipment reached stars of magnitude 13™:5 in the Hyades, with standard errors 
of o™-050 for the early-type stars and 0o™-058 for later types; in the Praesepe 
region the limiting magnitude was 13™-0, and the standard errors 0™-038 and 
0™-045 for early and late stars respectively. ‘The dispersion used was about 
200A./mm. and the exposure time for the faintest star was two hours. In 
comparing Ramberg’s standard errors with those in the present work, one must 
recognize that comparisons between stars well separated in the sky may be 
expected to exhibit considerably larger errors than comparisons in a restricted 
region photographed together, other things being equal. 

4. Eight of the stars in Table I have been observed by Mr Hogg and Miss 
Mason with a photoelectric cell placed behind an aperture in the focal plane 
of the slitless spectrograph on the 30-inch reflector. This aperture passes 100 A. 
of the stellar spectrum, centred on 45504. Comparisons are shown in Table II 
between relative magnitudes observed with the photocell (preliminary values, 
subject to revision) and those observed with the coarse grating, the latter being 
shown both uncorrected for the wave-length difference between 4550A. and 
4480. and corrected for this difference by subtracting 0-038¢, the gradient ¢ 
being taken from Table I. 
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TABLE I 
Observed Magnitudes 
Standard errors 
Mag. Mag. ,. No. of 
Star Spectrum (5980 A.) (4480 A.) (in — of . — plates 
Tucn Bg 6-4 
A Hydi Ks 
a Erid Bo 
n Circ Gs 
e Circ Ko 
y Tr Au Ao 
e Tr Au Ko 
B Tr Au Fo 
6 Tr Au 
¢ Tr Au 
y Apus 
B Apus 
a Tr Au 
« Apus 
5 Arae 
7 Pavo 
Octn 
e Sgtr 
Pavo 
Pavo 
A Pavo 
€ Pavo 
3 Pavo 
B Pavo 
Octn 
o Pavo 
y Pavo 
v Octn 
Octn 
« Tucn 
v Indi 
B Octn 
y! Octn 
y? Octn 


Comparison of Coarse Grating Magnitudes with Photocell 
Coarse Photocell Co 
Photoelectric ; A grating minus 
grating minus 
4550 A. corrected corrected 
4 =e to 4550 A. grating 
m ™m m 


m™m m 
a Erid Bo —1°47 —1°52 +0°07 +0:06 
5 Arae Bs +1°72 +1°93 
e Sgtr Ao 0°00 + 0°10 +o'10 
y Tr Au Ao * + 1°04 +0°03 + 1°04 +o°o1 
B Pavo AS +1°76 —0'05 +1°75 —0'04 
B Octn Fo +2°44 0°00 +2°43 
5 Pavo Gs +2°23 +2-28 —0°05 +2°23 0-00 
a Tr Au K2 +1°27 —o'1o +1°17 


Star Spectrum 


The root mean square of the “ uncorrected ’”’ residuals is +0-06,, and that 
of the “corrected”’ residuals is +0-04,. The agreement is very satisfactory 


HM 
5°32 —oro2 
4°76 +2:27 
0°60 
+1°53. 
4°68 +2:-01 
2°98 +0°02 
3°95 
3°93 +196 
4°83 
3°74 
4°00 
1°72 +2°47 
5°49 0°00: 
3°68 —o-18 
5°53. + 0°05 
3°94 +2°08 
5°84 +0°56. 
4°28 
4°02 
= 
517 +0°79 
4°84 
H 3°59 +1°81 
2°74 +2°54 
5°33 +1°33: ie 
4°28 +0°31 
4°99 +1°69 
5°52 
Tasxe II 
| 
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and shows that there is no difference in scale over 4 magnitudes and also (since 
a wide range of spectral type is concerned) that the blue grating magnitudes are 
very close to being monochromatic magnitudes, despite the low resolution. 

Eight of the stars have been measured in Dr Gascoigne’s relative gradient 
programme, also carried out on the 30-inch reflector. The relative gradient 
between two stars is defined by 


_ I. 4m,—Am, 
$= — I/A,)" 


Taste III 
Comparison of Gradients 

Diff. Diff. 
Star Spectrum 4, Gascoigne ¢, Coarse grating Gascoigne minus minus 

Coarse grating 0°08 
« Erid Bo —0'24 —0°35 +0°03 
A Pavo B2 —or'1o +0°05 —0'03 
6 Arae B, —o'l4 —o'18 +0°04 —0'04 
e Sgtr Ao +004 +0°07 —0'03 
y TrAu Ao +0'03 +o-o1 —0'07 
B Tr Au Fo +0744 +024 +o°16 
B Pavo As +0°37 +o°12 +0°04 
Pavo Gs +1743 + 1°37 —o-02 


The mean difference is (Gascoigne minus Coarse Grating) +0-08 indicating a 
difference in the adopted zeros of the scales. The differences between the relative 
gradients are shown in the last column, the root mean square being +0-07,. As 
the standard error of the Greenwich gradient determinations is + 0-03, this 
comparison is quite satisfactory. 

5. Fig. 2 shows gradients plotted against spectral type, and Fig. 3 shows the 
difference between the magnitudes interpolated to 5500. and the Harvard 
Mimeograph magnitudes plotted against gradients. There is a very marked 
dependence of the difference Stromlo minus HM on spectral type. This 
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Fic. 2.—Gradient plotted against spectral type. 
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Fic. 3.—Comparison with Harvard Mimeograph magnitudes. The magnitudes interpolated 
to 5500 A. are shown plotted against the gradient 4, various symbols representing the spectral type- 


MKGFAB 


difference is actually decreased if the earlier Harvard magnitudes HR are 
substituted for the more recent ones. In the case of these bright, far southern 
stars, the HR magnitudes are simple means of Bailey’s visual determinations 
with a 4-inch photometer (Harvard Annals, 34 and 46). The corrections applied 
to the magnitudes in HR to bring them to HM are shown in Table V of Miss 
Payne’s Photovisual Magnitudes in Harvard Standard Region (Harvard Annals, 89). 
These corrections were derived by comparing Bailey’s visual magnitudes for 
g8 stars between 6™ and 9™ with standard photovisual stars. Miss Payne’s table 
shows that Pickering’s colour correction with the 12-inch photometer was 
strongly dependent on magnitude between 7™ and 9g”, and it is suggested that the 
adopted corrections to Bailey’s observations of fainter stars are not applicable 
to brighter stars. This conclusion is supported by Cousins *, whose photovisual 
magnitudes of stars brighter than 6™-o show a strong colour difference with HM. 

6. The magnitudes and gradients, determined with a low resolving power 
such as this, might be called coarse monochromatic magnitudes and coarse 
gradients, to distinguish them from quantities observed with instruments having 
higher resolution. It is however undesirable to introduce new technical terms 
when there is small justification, and moreover it is not certain where one should 


* A. W. J. Cousins, M.N., 103, 154, 1943- 
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draw the line between a coarse monochromatic magnitude (or oligochromatic 
magnitude) and a “‘ monochromatic ’’ magnitude. If the spectra were free from 
absorption lines, the simple and obvious definition of a monochromatic magnitude 
would be useful: but since stars do show absorption lines, for many purposes 
one does not wish to ascertain the light flux in a narrow range of frequencies which 
might fall within a line. The average value of this quantity over a range of 
frequencies large enough to sample the incidence of lines fairly may therefore 
be of far more value. We do not wish to labour this difficulty, but think it of 
sufficient weight to be considered when contemplating the introduction of 
qualifying adjectives to the term ‘‘ monochromatic magnitudes ”’. Regarding the 
coarse gradients, we feel that they are much closer to (accurate) gradients than 
they are to colour indices, because the spectral ranges considered are much more 
definite than those employed by instruments which take in a wide enough spectral 
range to be appreciably influenced by variations with wave-length of transmission 
in the Earth’s atmosphere and in glass, or in the reflectivity of silver. 

While, then, the authors expect the magnitudes and gradients determined with 
the low resolution they have used to show (small) differences from determinations 
made with higher resolving power, they have not thought it advisable to introduce 
special qualifying nomenclature. 

Acknowledgments.—The authors are particularly indebted to Mr J. H. 
Banham for the care with which he has supervised the construction of the 
apparatus, and to Miss Janet Lamb for measuring many of the plates. They also 
acknowledge with pleasure their indebtedness to Mr Hogg, Dr Gascoigne, 
Mr Stibbs and Miss Mason who have helped in a number of ways. 
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ON THE INTERPRETATION OF MEASURED SOLAR 
WAVE-LENGTHS 


Lyman Spitzer, Jr. 


(Received 1950 January 13) 


Summary 


Reasons are given for expecting that collisions with H and He atoms 
produce a shift of most Fraunhofer lines to the violet rather than to the red. 
When allowance is made for this effect, and for the many uncertainties in 
collisional broadening theory, the interferometric measures of solar wave- 
lengths seem consistent with a complete absence of systematic Doppler shifts 
on the solar surface. 


In a recent thorough-going and important study, Miss Adam (1) has measured 
with an interferometer the precise wave-lengths of 14 lines in the solar spectrum 
at 14 distances from the centre of the disk. For cos @ equal to 0-181, com- 
parison between solar and laboratory wave-lengths shows an average red shift of 
0-0129A., the exact value predicted by general relativity theory. With increasing 
distance from the limb the red shift was found to diminish, and at the centre of 
the disk amounted to less than half of the relativity shift. 

In the theoretical interpretation which Miss Adam develops for these 
observations, a substantial correction is made for the shifts produced by collisions 
between absorbing atoms and atomic H. The line shift, assumed to be toward 
the red in all cases, is set equal to 0-36y, where y is the collisional damping 
constant and the factor 0-36 is from a theory by Lindholm (2). Values of y at 
different distances from the limb are taken from the observational study of strong 
Fe lines by ten Bruggencate and Houtgast (3). Since y is greatest at the centre 
of the disk, where we see down into deeper layers in which collisions are more 
frequent, the theoretical collisional red shift is greatest near the centre of the Sun. 
Since the observations show a red shift which is actually least at the centre, 
appreciable upwards velocities, at least 0-5 km./sec., seem indicated. However, 
such velocities are too great to be reconciled with existing theories of the solar 
atmosphere. 

The theory of collisional perturbations on absorption lines plays a major 
; part in this analysis. The following brief summary of relevant information on 
. this subject indicates that certain modifications seem appropriate in the theory 
: used by Miss Adam. 

7 The shift of absorption lines in the presence of atoms of a foreign gas has 
i been known for some time. Observations and theories in this field have been 
i summarized by Weisskopf (4) and by Margenau and Watson (5). This latter 
' survey also presents the previous theory of Margenau, which accounts for the 

observed fact that the shift of a broadened line averages a little less than one-half 

of the half-width. The subsequent contribution by Lindholm (2) lies not in 
discovering this proportion between shift and half-width, but rather in presenting 
a simplified form of collisional broadening theory, in which the ratio of collisional 
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shift to collisional half-width may be computed for any assumed law of interaction 
between the colliding atoms. This same theory has also been presented 
independently by Foley (6). When the perturbing force varies as the inverse 
sixth power of the separation, as may be expected for interaction at great distances 
between H and other atoms, this theoretical ratio becomes 0-36. 

It should be emphasized that this theoretical result is only an approximation. 
In the first place the theory, while a considerable improvement over previous 
theories, is still drastically over-simplified. A discussion of the theoretical 
proble: is involved is given in Unséld’s survey article (7) on collisional effects in 
spectral lines. Moreover, the observations, which have been summarized by 
Foley (6), show a very large scattering of this ratio about the theoretical value, 
with values ranging from 0-16 to 0-77. While the average ratio is 0-38, in close 
agreement with the Lindholm-Foley theory, any use of a constant ratio involves 
the possibility of very considerable error. 

The direction of shift, obviously a matter of great importance, will depend 
on the difference in perturbation between the upper and the lower levels of the 
transition producing the line in question. If the force between the two atoms is 
attractive, the energy levels will be depressed. ‘The usual van der Waals forces 
will depress the upper level more than the lower one, and the line will be shifted 
to the red. This corresponds with the observed fact that almost all resonance 
lines are shifted to the red by the presence of a foreign gas. 

For levels of high excitation the van der Waals theory is not applicable, since 
the radius of the electron orbit becomes so large that a perturbing H atom passes 


inside this orbit during a typical collision. The theory for a highly excited level . 


(n equal to 15 or more) has been considered by Fermi (8) and Reinsberg (9). The 
computations show that the sign of the dominant term in the perturbation energy 
depends on the phase shift 5 of the wave function for a free electron when such 
an electron is scattered by the perturbing atom. If tand is negative, the energy 
level is raised, with a violet shift resulting. Observations show that the higher 
members of the principal series in the alkali metals are shifted to the violet by 
He and Ne, and to the red by A, Kr and Xe. The discussion by Reinsberg (9) 
indicates that these results are in good quantitative agreement with the theory. 

While there are no available observations on collisional shifts of this type 
produced by H, the theory may be applied in this case also. Wave functions for 
a free electron in the field of an H atom have been computed by Morse and 
Allis (x0) ; electron exchange was taken into account, but separable wave 
functions were used. A simplified form of non-separable wave function, which 
neglects electron exchange, has been considered by Huang (11). Both of these 
results show a phase shift decreasing monotonically from 180° at zero energy 
to a low value at high kinetic energies, a result similar to that found for He. 
Thus tan3 is negative for electrons of low energy. It follows that collisions with 
H atoms will produce shifts to the violet in any absorption line to a highly excited 
state, near the ionization limit. This result does not seem to have been realized 
previously. 

For lines arising from transitions between levels of intermediate excitation 
there is at present no adequate theory available. Observations on the spectra 
of alkali metals indicate that where violet shifts appear at all they appear in all 
lines except the resonance line. Specifically, when He is the perturbing gas, the 
second and all higher members of the 6s—np series of Cs (12) and the second mernber 
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of the 5s—-np series of Rb(13) are shifted tothe violet ; for K (14)the measures extend 
down only to the third line of the 4s—np series, but the violet shift of this line 
exceeds the shift near the series limit. Thus violet shifts of most absorption lines 
seem probable when an atom is perturbed by H or He. 

While collisional violet shifts seem likely in the Sun, no precise theoretical 
prediction of such shifts is presently possible. One may nevertheless ask whether 
Miss Adam’s measured shifts are consistent with the belief that no systematic 
Doppler shifts are present, and that the observed values of A, the red shift 
measured in units of 0-001 A., are the combined effect of collisional shifts and the 
Einstein gravitational shift. 

To answer this question we first repeat Miss Adam’s analysis of her data, 
changing the sign of the collisional shift but retaining all her other assumptions. 
The apparent Doppler shifts resulting are given in Table I below, which follows 
closely the form of Miss Adam’s Table V. It is evident that the resultant 
velocities are now only about one-third of those previously found. 


TaBLe | 
Mean Line Shifts and Radial Velocities 


Collisional Residual Velocity shift Velocity vg 
and p.e. km./sec. 
0°639 —4°3+0°6 —O'211 +0°029 
0°468 —~4'1 +0°6 —0'201 +0°031 
0°398 ~3:1+06 —0'152+0°032 
0-282 —1'7+0°7 — +0°032 


cos 0 


Three additional modifications must also be considered. Firstly, there 
may be a systematic error either in the solar wave-lengths or in the laboratory 
values. Since the observed value of A at cos@=0-181 is just equal to the value 
predicted from general relativity, it is evident that this value of A must be in error 
if collisional shift is assumed present without any Doppler shift. Secondly, 
the damping constant for the lines investigated need not be the same as for the 
lines studied by ten Bruggencate and Houtgast (3), since considerable scatter 
of damping constant may be expected from one line to another because of 
variation with atomic properties. Thirdly, the ratio between collisional shift 
and damping constant may differ markedly from the theoretical value. With 
so many uncertainties in the theoretical interpretation, the resultant Doppler 
velocities can evidently be reduced to a low value. 

On the basis of the evidence available one may conclude that Miss Adam’s 
observations are consistent with the belief that motions in the solar atmosphere 
(excluding rotation and small-scale phenomena near sunspots) are so slight 
that the measured shifts are simply the sum of the relativity red shift and a 
collisional violet shift. On the other hand, these shifts are not inconsistent with 
moderate Doppler shifts of some 0-2km./sec. As suggested by Richardson and 
Schwarzschild (15), effective shifts of this order could be produced by rapid 
convective currents with a scale much smaller than the observable granules. 
The assumption that both Doppler shifts and collisional shifts are present is 
consistent with the extensive observational material obtained by Adams (16) 
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on the red shift of Fraunhofer lines at the Sun’s limb relative to the centre. The 
relatively low shifts found for lines formed high in the Sun’s atmosphere, where 
large turbulent velocities might be expected, would seem somewhat unfavourable 
to the hypothesis that Doppler shifts alone could account for the centre-limb 
variations. 


Princeton University Observatory : 
1949 December 23. 
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‘THE CAPE AND RADCLIFFE OBSERVATIONS OF THE 
ECLIPSE COMET (19481) 


F. v. B. Lourens and T. W. Russo 


(Communicated by H.M. Astronomer at the Cape of Good Hope) 


(Received 1950 January 2) 


Summary 


Positions are given for the Eclipse Comet 1948 | derived from plates taken 
at the Cape Observatory with various telescopes between 1948 November 8 and 
1949 March 2, and plates taken with the 74-inch reflector of the Radcliffe 
Observatory, Pretoria, between 1948 December 2 and 1949 April 2. All the 
plates were measured and reduced at the Cape Observatory. 


Table I gives the observed position of Comet 19481 on 41 occasions between 
1948 November 8 and 1949 April 2. The first three columns give the universal 
time of the observations and the measured positions reduced to the FK3 system 
and referred to the equinox at the beginning of the year. The three columns 
headed A, B and H give the number of exposures measured, the number of 
reference stars used and the hour angle of the comet, while the sixth and final 
columns, respectively, give the weight and the estimated magnitude, assigned to 
the observed position, unit weight corresponding to a probable error in position 
of +0"-3. The positions given still require corrections for :— 


(1) aberration due to the Earth’s motion, 


(2) the:time taken by light in passing from the comet to the Earth, 
(3) parallax. 


The correction for parallax can be deduced from the quantities given in the 
last two columns but one of the table by dividing them by the geocentric distance 
of the comet. The tabulated quantities were computed from the formulae 


A. p, = +0*- 4873 sin H sec 3, 
A. ps = — 4"°885 cosd — 7”-309 sind cos H 


for the Cape observations and from 


A.p,= +0*5287 sin H 
A. ps = — 3":806 — 7"-931 cos H 
for the Radcliffe observations. 

The observations at the Cape were made with three instruments: the 5-inch 
Wide-Angle Astrographic Camera (scale 1 mm. =102"-52), the 13-inch Astro- 
graphic Refractor (scale 1 mm.=59"-81) and the 24-inch Victoria Refractor 
(scale I mm. =30"-62). From November 8 to December 13 the comet was 
photographed with both the first two of these instruments and one plate 
(containing several exposures) from each instrument was measured for each night, 
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TABLE | 


Cape Observations: (Long. —1® 13™ 558, Lat. —33° 56’ 03”) 


Date (U.T.) 


R.A. 


Dec. 


A B Wt. 


H 


A. da 


1948 


8°11054 
10104 
10°09404 
1309287 
14°10464 
15°08394 
17°09345 
19°07219 
22°07665 
24°09747 
26°08575 
27°06825 
28-06398 
2704484 
8-05688 
13°07266 
14°04005 
15°04456 
21°87236 
29°88343 
30°93265 
1949 
5°91450 
9°90193 
10°89003 
23°84553 
31°81252 
3°83591 
482225 
582338 
780940 
15°79840 
1680959 
17°84287 
2°79542 


Nov. 


hm 
13 25 21°1 
13 17 22°4 
13 9 23°1 
12 45 13°84 
1237 4°44 
1229 9°83 
12 12 52°59 
11 56 43°80 
II 31 56-90 
15 4°36 
10 58 19°14 
10 49 59°66 
10 41 31°80 
10 7 36°56 
9 17 25°89 
8 38 27°73 
8 31 23°92 
8 24 14°33 
7 40 46°42 
7 1 33°54 
6 57 18-31 


6 36 33°92 
6 25 36°57 
6 23 12°52 
6 0 57°02 
5 53 45°09 
5 51 58°40 
5 51 29°68 
551 2°96 
5 50 18°33 
5 48 52°11 
5 48 50°59 
5 48 51°00 
5 51 19°91 


1948-0 


—21 18 35°0 
—22 6 44°0 
—22 54 16°3 
—25 12 43°0 
—25 57 22°4 
—26 39 29°8 
—28 
—29 18 10°7 
1 5§2°0 
—32 23671 
—32 54 30°2 
—33 17 29 
—33 37 44°4 
—34 37 15"! 
— 34 53 14°2 
—34 13572 
—33 45 40°4 
—33 27 23°2 
— 3° 44 39°5 
—26 44 45°6 
—26 11 54°0 


1949°0 


—23 5 308 
—21 §20°1 
—20 36 20°8 
—14 §2 25°9 
—11 §5 21°9 
—10 54 37°3 
—10 35 34°0 
—10 16 31°6 
— 939 51'4 
— 725 32°3 
—7 9558 
— 6 54 163 
393 


AL AAU AWS 


S 


h m 
—6 23°6 
—6 25°4 
—6 23°6 
“SOs 
20°2 


| 4 


—4 23°9 
—3 29°1 
—3 21°4 
—3 34°4 
—3 
—3 671 
—1 
13°4 
49°4 
—o 
—3 29°4 
—2 29 
43°4 


—025'1 
—o 16°6 
27°4 
—o 18:2 
—027°1 
+0 20°3 
+o 5°0 
+o 
o's 
+o 16°5 
+0 36°7 
+1 28-7 
+1 89 


Radcliffe Observations : (Long. —1% 52™ 55%, Lat. —25° 47’ 18”) 


1948 
299954 

1949 
22°79889 
3°76212 
15°75890 
23°75984 
30°80538 


Apr. 2°80052 


9 59 30°1 


5 49 18°47 
5 51 40°32 
5 57 12°76 
6 2 137 
6 6 49°60 
6 9 0°27 


1948-0 


—34 45 47°7 


1949°0 


— 5432271 
— 35° 47°5 

I 47 33°6 
— O40 11'S 
+ 010 36°8 
+ 0 29 584 


—3 20°3 


+1 23°4 
+1 3°5 
+1 40°6 
+2 86 
+3 37°72 
+3 39°8 


—0°494 


+0189 
+0°145 
+0°225 
+0281 
+0°429 
+0°433 


— 3°05 
—3°28 
—3°58 
—3°73 
—3°82 


except on November 9 when only the Astrographic Refractor plate was measured 
and on November 13 when two plates from each instrument were measured. 
For the first three observations the comet was very bright and only 7° above the 


—o'520 —4°82] 1°8 
—0o'523 —4°83] 2°0 
—o'538 —4':27| 2°7 
—0'534 —3°84] 28 
| —5 38:2 | —0°543  —4:05| 3°0 
o-3 |  —3°43| 3°4 
1'0 —0'544 —3°44] 3°8 
—o'519 —2°65/| 4°3 
1'2 —0'455 —1'°77| 4°6 
I'l —0°447 —1'57} 
—0o'469 —1°70] 5°2 
Dec. 12 —o'430 ~—1'16] 6:0 
—o'239 —o'18] 6°9 
I'l —0'034 +0°03|] 7°6 
—o'081 —o-08| 7:8 
—0o'449 8-7 
—o'278 —1'53] 9°7 
1-0 —o-'102 —1'21| 98 
Jan. 10 —o'058 —1°64| 10°4 
10 —0-062 10°9 
08 —or0g0 | 1271 
—o'059 | 12°8 j 
Feb. +0°044 —3°42] 13°'0 
| +0024 —3°50| 13°2 
+0°035 —3°90| 13°8 
+0078 =—3'95 | 13°9 
+0185 —4°04]| 14°0 
Mar. +0145 —4°38] 14°5 
Feb. a 4 eg 14°2 
Mar. 3 14°5 
| 3 03 | | 
— 24 05 | —3°85 | 16-0 | 
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horizon. Accurate measurements of its position were therefore very difficult, 
especially in right ascension, as the direction of the tail was mainly east to west. 

From December 14 onwards the observations were made with the Victoria 
Telescope. Up to December 14 the telescope was guided on the comet and the 
stars allowed to trail; after December 14 a star was used for guiding and the comet 
allowed to trail. For the last four Cape observations neutral filters were used to 
reduce the brightness of the comparison stars. These filters were prepared by 
fogging photographic plates and then removing the film from a small area in the 
centre through which the comet was exposed. One such filter reduced the apparent 
brightness of the comparison stars by 2-8 magnitudes and two filters by 5:5 
magnitudes. An exposure of from 14 to 30 minutes was required to record the 
comet on Ilford Astronomical Zenith Plates which had been sensitized with 
mercury vapour. 

The observations at Pretoria were made by Dr A. D. Thackeray with the 
74-inch Radcliffe Reflector (scale 1 mm. = 22"-50) and were reduced at the Cape 
Observatory with the aid of faint reference stars having photographic magnitudes 
between 11™ and 13™ whose positions had been determined from plates specially 
taken for this purpose with the Cape Refractors. The Radcliffe exposures were 
also guided on the stars. ‘The images of the comet are rather faint and difficult 
to measure. 

When two plates of the comet were available for the same night, one was 
reduced by the usual dependence method using three stars, the other by standard 
coordinate computation and the use of four quadrantal stars, the reference stars 
for the two methods being selected as differently as possible. When one plate 
only was available for a night, this plate was measured twice, once for each 
method. On November 17, however, only one reference star near the comet was 
used. 

The positions of the reference stars were taken from the Yale Photographic 
Zone Catalogues when the comet was north of — 30° declination and from 
unpublished Cape Photographic Zone Catalogues when the comet was south of 
— 30°. 

The epoch of observation of all the reference stars is comparatively recent, 
i but as proper motions were available their positions were reduced to the epoch — 
of observation of the comet. The positions given in the Cape Catalogues are 
already on the FK3 system. Those taken from the Yale Catalogues were reduced 
to the FK3 system by the small corrections shown in Table II. These quantities 


: ; were obtained from the differences ‘‘GC minus Yale” given in the introduction 
Taste II 
R.A. Yale Zones FK3—Yale 
13 —20 to —22 —o-008 +0o°21 
| 13 —22 to —27 —o'006 ++. 0°17 
i 12 —27 to —30 —0'002 0°23 
7 —22 to —27 — 
6} —20 to —22 +o010 +019 
6 —14 to —18 —0o'17 
6 —10to —14 —o'004 —0'28 
6 — 6to —10 —o'009 0°30 
6 — 2to— 6 —o019 +0°39 
6 —2to+ 1 —0:086 
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of the various Yale Catalogues together with the differences “‘ FK3 minus GC” 
as given in the G.C. and also by Kopff in A.N., 269, 3, 1939. 

The estimated magnitudes are based on a graphical smoothing of 40 direct 
estimates (to the nearest half magnitude) made between 1948 November 8 and 
1949 March 30. An investigation of the empirical relation between the observed 
magnitude (m) and distances (r, A) of the comet yielded the following result :— 

m = 6-97 + 5°76 x 2-5 logr + 5 log A, 
so that the magnitude of the comet when at unit distance from both Sun and Earth 
was 7:0, and varied approximately as the sixth power of its inverse distance from 
the Sun. The exponent given, 5-76, is an average over the range from r =0-5 to 


_r=3°0. Its actual value varied from 6-3 at r=0-5 to 5-4 at r=3-0. 


Royal Observatory, 
Cape of Good Hope : 
1949 May 14. 
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THE EXPANSION OF AN INTERSTELLAR GAS-CLOUD 
INTO A VACUUM* 


G. C. McVittie 


(Received 1950 January 30) 


Summary 


Riemann’s work is summarized for the case of the adiabatic one-dimen- 
sional motion of an interstellar gas-cloud expanding intoa vacuum. Boundary 
conditions are found for: (A) a homogeneous (isothermal) gas-cloud; 
(B) a non-homogeneous cloud having a temperature-gradient perpendicular 
to the face advancing into the vacuum. Due to the asymmetry of the 
boundary conditions, it is not possible to prove that one of the Riemann 
variables, s, is constant in case (A). General solution of the problem for case 
(A) under assumption that s is constant and relation to Burgers’ particular 
solution. Demonstration that atoms of the gas at the advancing face must 
be given an initial impulsive velocity in case (A). General solution for case 
(B), in which neither Riemann variable is constant. Advancing face shown to 
be permanently the seat of “ shock” conditions. Method of fitting the 
solution to any law of temperature-variation in the cloud given a priori. 
Particular example worked out in detail which suggests that Burgers’ high 
speed of advance of the face into the vacuum is attainable only after an 
infinite time. 


We consider a gas-cloud with an infinite plane face, the cloud extending to- 
infinity to the left of the face whilst to the right there isa vacuum. The gas-cloud 
consists of monatomic hydrogen so that the ratio of specific heats, k, is §. The 
equation of state is p=(R/m)pT and the temperature of the gas is about 10,000 
deg. K. ‘The motion is one-dimensional and the x-axis is taken perpendicular 
to the face. Adiabatic conditions hold throughout the motion so that p=«xp* 
and there is no viscosity, heat-conduction, radiation pressure, etc. A solution 
of the problem under these conditions has been given for a homogeneous gas-cloud 


by Burgers +, who arrives at the somewhat surprising conclusion that the face of — 


the cloud advances into the vacuum with three times the speed of sound in the 
undisturbed gas. This is about 35km./sec., which we may compare with the 
peculiar velocities of some 15 to 20km./sec. of the clouds themselves. The 
object of the present investigation is to show that the problem is much more 
complex than Burgers’ treatment would suggest, and to indicate alternative 
methods of tackling it. 

1. General Theory.—The fundamental equations governing one-dimensional 
supersonic flow may not be generally familiar to astronomers and we therefore 
summarize in this section Riemann’s{ basic investigation of this problem. 

* A summary of this paper was read at a Joint Colloquium on Problems of Motion of Gaseous 


Masses of Cosmical Dimensions, Paris, 1949 August. 


+ J. M. Burgers, K. Ned. Akad. v. Wet., 49, 2, Section 2, 1946. The physical data on 
interstellar gas-clouds are summarized in Section 1 of this paper. See also J. H. Oort, M.N., 
106, 159, 1946, and W. S. Adams, P.A.S.P., 60, 174, 1948. 

t B. Riemann, Ocuvres Mathematiques, Paris, p. 177, 1898. 
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‘The equations of motion and of continuity in the usual notation of hydrodynamics 


-are 


pot pdx dx 
If the adiabatic condition holds, we can replace p by the velocity of sound c, where 


oO. 


=kp/p =kxp* =k(R/m)T, and the equations become 


Ou = 2e Oc 
+ =0, (1.1) 


2 Oc 2u dc Ou j 
+ =0. (1.2) 


‘The Riemann variables r,s, where 
(1.3) 


(1.4) 


transform (1.1), (1.2) into 


(1.5) 


(1.6) 
But by (1.3), (1.4) 


u+c=ar+ Bs, u—c=—fr—as, (1.7) 
where 


a=4(k+1) a+B=k—1 
or 
B=Kk-3) a—f=2. 
Hence (1.5), (1.6) become 


(1.8) 


(1.9) 


(1.10) 


we can transform to r, s as independent variables. In fact, since 


dx, 


* B. Riemann, /oc. cit. 


| 
Ou Ou 10p 
at ~~ pax’ 
1 Op u Op | du 
or 
2c 
or or 
|| 
ar ar 
: + (ar + Bs) = = 0, | | 
ds ds 
az (Br+as) = 0. | | 
It has been shown by Riemann * that, when | | 
or Os 
7” | 
‘ 
= 
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we have by (1.9), (1.10), 


Bs) dt + dx} = Bs)t) + ta dr + Bds)}, 


and similarly 


Hence if we write 


X =x—(ar+ Bs)t, Y=x+(Pr+as)t, (1.11) 


and regard X, Y as functions of r,s, we have 


Or { (0X ox 
dr = + at) dr + (5 + 


as (/aY oY 
ds = pt) dr + at) ds}. 


But r,s are independent variables, so that the last two equations give 


(1.12) 


Since, by (1.12), we have 0X/0s —(—dY/0r) =o, we can define the total derivative 


of a function w(r,s) by a 
dw =X dr—Y ds (1.14) . 


or, using suffixes to denote the partial derivatives of w, 
w,= —Y=—x—(Br+as)t, 
w,=X =x—(ar+ Bs)t. 


(1.15) 


The equations (1.12) and (1.15) yield 


(1.16) 


Hence 

n 

Wr, + Ws) =O, (1.17) 

where 
n= a+B (1.18) 


For an arbitrary gas, » may have a non-integral value, but for a monatomic 
gas n=I, a result with important consequences, as we shall see in Section 5. 
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The analysis leading up to (1.17) presupposes that both 07/@x and 0s/0x are 
not zero. If one or other of them is zero, it is best not to introduce the function w 
but to solve the equations (1.9) and (1.10) directly. We note also that the analysis 
breaks down when either or both of @r/0x, ds/0x is infinite. Regions in which this 
happens we shall describe as being the seat of “shock” conditions: in them the 
equation p=xp* no longer holds good. 

2. Boundary Conditions.—Suppose that initially the face of the gas was at 
(B,);, where x =(x,);, and that ata later time it is at B, and that B, is then advancing 
to the right (see Fig. 1) with speed v,. The cloud then lies to the left of B, whilst 


B. (= 8: i 
GAS | IN VACUUM 


| 


x= 


Fic. 1. 


to the right of this plane there is a vacuum. Following Burgers, we define B, 
as a plane on which the density of the gas is zero, i.e. by (1.3) 


+5, =0. (2.1) 


The speed of advance of B, is v; =(dx/dt),, where, by total differentiation of (2.1), . 


ar 
dx * dx), \at), 


Or as 
(+5), * 


Using (1.5), (1.6) with u=u,, c=0, we obtain 


(2.2). 


whether one or other of (0r/0x),, (s/0x),.is zero or not. ‘Thus the boundary B, 
advances with the speed of the gas at B,, as it should do. 

Strictly speaking (2.1) is the only boundary condition, for we do not 
know how the expansion started, nor do we know how the successive layers of gas 
are set in motion as time progresses. However, it is reasonable to suppose that all 
the gas does not start moving at once and that there is an inner boundary, B,, to 
the moving part of the gas. To the left of B, the gas is at rest and the pressure 
and density are continuous from one side of B, to the other. Initially B, is at 
(B,);, where x =(x);, but we do not for the moment commit ourselves to the state- 
ment that (B,);=(B,);. The conditions at B, are ug =0, Cc, = — Cy or by (1.3), (1.4), 


(2.3) 


9+ — 2€9/(R—1). (2.4) 


SAS 
aT 
RES] 4 

X=X2 

i 
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Case (A). Homogeneous (isothermal) cloud.—In such a cloud ¢g is a constant. 
Let B, be at x, at time ¢, and at x, — dx at time t,+ dt. Then the speed of advance 
of B, is v,=(—dx/dt),, where, by total differentiation of (2.3), (2.4), 


(5 Ox Ox dt 


+ + 5) (=) 
ot Ot] 0. dt 


Using (1.5), (1.6) with u=0, c= —¢9, we have 


Ox Ox ax] 


whence we obtain the two conditions 


(2.7) 
Hence B, moves to the left with the speed of sound in the undisturbed gas. 
As to the solution of (1.9), (1.10) which will give the motion in the region from 
B, to B,, all we can say is this:—Let (B,),; be the initial position of B,. Then 


at each value of x lying between (x,);and x, and so s is constant for such values of x. 


But this gives us no clue as to the way s varies for values of x between (x); and x, 
even if initially (x.);=(x,);. No doubt we can assert that s is constant for the whole 
range x, to x, on the basis of some additional assumption, e. g. that s is an analytic 
function of x. At any rate the boundary conditions alone are insufficient for the 
purpose and this arises from their asymmetrical character. Had, for example, B, 
been advancing into the same gas at rest, it would have been easy to show that s 
was constant in the region from (B,); to B, also, but this is not the case in our 
problem. 

Case (B). Non-homogeneous cloud.—We shall now show that the conditions 
(2.6), (2.7) depend on the postulated homogeneity of the cloud. .Let us suppose 
that the cloud is non-homogeneous in such a way that cy=Cc (x) is a function of 
x given a priori and valid for that portion of the cloud which, at each instant, 
has not yet begun to move. We can interpret this condition physically by saying 
that there is a temperature-gradient in the cloud in a direction perpendicular to 
the face. We write c, for dcy/dx and differentiate (2.3), (2.4) totally again. We 


find 
Or as Or as —dx\ _ 
(5 si) (5 FF. 
ot ot 2 Ox 2 dt 2 —1\ dt 2° 
whence 
Or as Or as 
a= (32) ax Ox 
os| | “dr os 206 


ve ( 
ds 
Ox 2 
2 
a 
4 
a 
| 
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Thus we now have instead of (2.6), (2.7) 


Oxdx), k—1\° \dx dx] J,’ 


2 


We see at once that a solution of (1.9), (1.10) with 


Ox 
will not do for this case, whatever may happen in case (A). Moreover the non- 
linearity of the equations does not permit us to argue a priori that the motions will 
possess any feature in common. 
3- Boundary conditions in terms of derivatives of w.—In cases where 
Or Os 
Ox Ox 
it is convenient to express the boundary conditions in terms of the derivatives of w. 
At B,.—Let B, be at x, at time t,. Then since r,+s,=0 at B, we can write 


(1.15) as 
(w,), = 


(w,), 


#0 


(3-1) 
Hence at B, 

(w,+W,) =0. (3-2) 
‘This condition being fulfilled, the speed of advance of B, is by (2.2), (2.1), (1.4), 
@,=u,=2r,. Hence (3.1), (1.8) yield 

(w,)) hy, (3-3) 


which is the equation for 7, in terms of x,, 4. 
At B,.—Here r, =s, and so (1.15) now becomes 


— ¥2—(%+ B)rote, | 
=X_—(%+ B)roto. 


‘Thus there is no condition analogous to (3.2) at B, except at the peculiar instant 
t,=0, which can always be avoided by changing the origin of time. Assuming 
that this has been done, we may solve (3.4), to give 


(3-4) 


3-5 

2 

But = — and «+8 =k—1 and therefore the last two equations 

give 


W, + W, 
17 


229 
ar, 

Ox Ox 
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Again, (1.13), (1.15) and (1.16) yield 


or ds 


Hence in case (B) of the preceding section the boundary condition (2.8) becomes 
2(k—1)= Wry) (3.8) 
provided however that 
Lastly the formula (2.9) for the rate of advance of B, becomes 


Ws, — 


Weg Wy, — 
4. The solution s = constant for case (A).—We now enquire how far the particular 
solution of (1.9), (1.10) for which s = constant will fit this case. We can solve the 
equations for any value of k and therefore we shall not introduce the value k = $ in 
our formulae. The equation (1.10) is satisfied identically whilst the constant 
value of s is determined by the boundary conditions at B, (equations (2.3) and 
(2.4)) to be 
(4-1) 


Combining this with the boundary conditions at B, (equations (2.1) and (2.2)) 
and using also (1.4), we obtain 


=27,= — 28, =2¢)/(k—1) (4.2) 


for the rate of advance of B,. ‘The equation (1.9) is now 


and its solution is 
(4.3) 
where f(r) is an arbitrary function of r. Given this function, the equations 


(1.3), (1.4) and (4.1) determine u,c in terms of x,t by the simultaneous solution 
of the equations 


c 


The rate of advance of B, has already been found to be —¢ (equation (2.7)). 
Lastly we obtain by partial differentiation of (4.3), with respect to x, 


(4-4) 


= (x =(at+f'(r))?. (4.5) 


f 
20 
| 
0 Bey or 
t —T x 
i 3 4 
2c 
Ke 
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It is easily verified that Burgers’ solution* of the problem is obtained by 
putting f= in (4.3), (4.4), his values for u,c being 


( 


But his results, that the face B, advances to the right with speed 2cy/(k— 1) whilst 
B, advances to the left with speed cy, do not depend on these particular expressions 
for u,c: they are valid for any u,c which satisfy (4.4). Our formulae differ from 
his in the values assigned to (x,,¢,) and (x,,¢,) which Burgers denotes by (x’,t’) 
and (x”,t”), respectively. Using (4.4) with c=0, u=u, =2¢9/(k—1), we find 


2 
x, ty =f (4-7) 


whilst (4.4) with u=0, c= —¢g gives 


4): (48) 


which replace Burgers’ (x/t)' =2c¢)/(k—1) and (x/t)” = —¢ 9, respectively. 
We now consider the initial instant, ¢;, of the expansion and we denote all initial 
values by the suffix 7. Equations (4.2), (2.3), (2.4) give 


(1): =(— = = (R- 1), (4-9) 
(72); —Co(R-1) and (uy);=0. (4-10) 
The obvious way—though probably not the only conceivable one—of defining 


the initial instant of the expansion is by the coincidence of (B,); and (B,); in Fig. 1. 
This means that (x,);=(x,),; for t;=(t,);=(t2);, whence from (4.7) and (4.8) 


Frat (G25) + (- (4.12) 
Also by (4.5), (4.11) we have 


((&),), -1(- }] (4.13) 


Now consider Burgers’ solution in which f= 0. By (4.11), (4.12) this implies 
that 


(4.6) 


t;=0, 
and, by (4.13), (4.14), that 


Hence at the commencement of the motion we have two features worthy of comment, 
viz.: (i) the face of the cloud is the seat of “‘shock’’ conditions, and (ii) the initial 


* J. M. Burgers, loc. cit. 
17* 


Bente. 
| | 
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outward speed of the face is 2¢9(k—1). As Burgers himself remarks, such an 
initial state is somewhat artificial. ‘The condition (i) survives even if we use the 
solution (4.3) instead of Burgers’. For, since f is arbitrary, we may limit it by 
choosing ¢;=0 in (4.11) and this then means that fis an even function ofr. Thus 
f’is an odd function of 7 and therefore the two values of dr/0x given by (4.13), 
(4.14) respectively are either different, which is impossible, or singular, which 
again makes the face of the cloud the seat of “‘shock”’ conditions at the initial 
instant. It is perhaps not relevant in a hydrodynamical treatment to introduce 
a consideration drawn from the kinetic theory of gases. However, we may remark 
that, if V is the root mean square velocity of the atoms of the gas when its 
temperature is 7, then V?=3(R/m)T, whilst the speed of sound is related to T 
by c?=k(R/m)T. ‘Thus when k= $ we have c?=5V?/9 which combined with 
(4.9) gives (u,);=3¢9=/5V. Hence the initial speed of the gas at the face is 
\/5 times the root mean square velocity of the atoms, whereas, in the absence of an 
impulsive force at the face, we should not expect the initial speed to exceed V. 
Our general conclusion is that this difficulty, as well as the high speed of advance 
of the face (3c, for a monatomic gas), is due to the assumption s = constant which, 
as we have seen, we are not compelled to make by the conditions of the problem. 

5. The general solution for case (B).—In this problem we have to use a solution 
of (1.9), (1.10) for which both dr/dx and 0s/0x are different from zero, which we 
find by solving (1.17). The integration of this equation for an arbitrary value 
of m requires the use of Riemann’s method *, but for our monatomic gas n=1 
and we can employ the simpler method of Laplace.t We write (1.17) in the 
equivalent form 


the integration of which reduces to the successive integration of the two equations 


(5 I) 

ow 
a (5-2) 


From (5.1) we obtain W= F(s)/(r+s), where F(s) is an arbitrary function of s, 
and then (5.2) gives the general integral of (1.17), viz. 


(5-3) 


where R(r) is an arbitrary function of r and S(s)=fF(s)ds. Writing S’ =dS/ds, 
R’ =4dR dr, etc., the various partial derivatives of w are 


+ RO) _S(s)_ S(s)+ 
R" 2R 2S+R) Ss" 28’ 2(S+R) | 


S’'+R'  2(S+R) 


J 


* G. Darboux, Legons sur la Théorie Genérale des Surfaces, Paris, Vol. II, chap. IV, 1915. 
1 G. Darboux, loc. cit., chap. II. 
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The equations for r,s (and so for u,c) in terms of x,t are, by (1.15), (5.4); 


S’  S+R 

—x—(Br+as)t= 
S+R 


whilst by (3.7) with «/B =(k+1)/(k—3) = —2, we have 


=(r+s)[(r+s)?R’ —2(r+s)(S’+2R’)+0(S+ 
9s (5.6) 
—2(r+s)(2S’ + R’)+6(S+ 
We now apply the boundary conditions. At B, we have r,;=—s, and by 
(3-2), (5-4) 
Ri(r)+S(s) 2(S+R) 
r+s (r+s)? 


(5-7) 


Suppose that R(r), S(s) are functions having Taylor expansions in the region 
between B, and B,. Consider a plane B} for which r+s=e, where « is small. 
Then B;—- B, as «+0 and r->r, and we can write (5.7) in the form 


e—>0 e 6 


For this limit to exist we must have 
R(r,)+ S(—1) =0, R'(r,)-— S(-—1) =0. 
But the first of these conditions implies the second, provided that 
S(s) = —R(-s), (5.8): 


and, of course, R’(—s) now means dR(—s)/d(—s)=—dR(—s)/ds. If now we 
calculate the values of 0r/0x, ds/dx at B,, using (5.6), (5.8) and the same limiting 
process, we find 


(5-9) 


(=), = Lime [ + Rr) + | 


Hence both (0r/0x), (ds/0x) have infinite values at B, which is therefore the seat. 
of ‘‘shock”’ conditions throughout the motion. At Bj the two equations (5.5); 
become, in virtue of (5.8), 


x—(2r—ae)t =}R’(r) + Ofe), 
x—(2r+ Be)t + 


| 
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Thus at B,, where « =o and r, = }u, = $v, *, we obtain, using also (1.8), the single 
equation 


which gives the speed of advance of B, in terms of x, t,. 


We next turn to B,, where rz =5_ = — 3c¢9(x2). It is simpler to work in terms of 
r, rather than of ¢o(x.) and we shall also use for brevity the notation 


R(r.)=R., = R’,, ete., 
R(-—1r.)=R_, R'(—12)=R‘ ete. 


Firstly, we can obtain the position of B, at any time by (3.5), (5-4), (5-8), 


R'.- 
TE 


I 
make 


(5.11) 


—r(R', + R_)+(R,—R_)}. 


We shall see in a moment that the first of these equations is equivalent to a 
relation between the functions R and cy. It is the second which then gives x, as 
a function of t, and so fixes the position of B, at each instant. Secondly, we can 
find the speed of B, by (3.10), (5.4), (5.8) and the fact that cy= — 3rg, viz. 

R’, — +R“) (5.13) 
3(R, -R.)—3r(R +R 
Thirdly, we can show as follows that the boundary condition (2.8) in its form 
(3.8) enables us to find the functional form of R when we are given the function 

€o(X2). We have by (3.8), (5.4), 


475 dx, 


= 


Using (5.11) this may be written 


Id 
4 dx, dry dr ' 


|R(r.)+ » |. 


which on integration yields the first of the equations (5.12), viz. 


But now r, = — 3¢y(x2) is a(known) function of x,, so that the last equation becomes 


d 
de, + 


4% 


* I am indebted to Professor E. T. Copson, University College, Dundee, for a proof of this 
statement. We have, if Bj — B,, that 


Since v,=dx,/dt, and u,=27r,, we obtain, on differentiating with respect to t,, 
vy =u, + f2t, + $R"(r,) dr /dty). 
But by (1.16), (5.4) and (1.8) with k=§, we have 
S+R S’+R’ 
(r+s)8 (r+s)?’ 


and if in this we put r=r,, s= —r, + €, use (5.8) and allow ¢€ —> 0, we have 2t,+4R”(r,;)=0, so that 
Uy. 


2 
(5.12) 

2 
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which on integration of the left-hand side by parts yields 
R(r2) + R( — 12) =2{C + — dxg}, (5-14) 


where C is a constant of integration. If cj! denotes the function inverse to Co 
we have x,=cj'(—%r,) and therefore we may express the right-hand side of 
(5.14), after performing the integration in the last term, as an explicit function of 
vz. We can then, in principle at least, select the function R which will satisfy 
(5.14) and thus determine the solution of (1.17) appropriate to the given 
temperature-distribution in the undisturbed cloud. 

We shall not discuss the initial instant of the expansion in this general case 
as we can illustrate the question more clearly by means of a specific example to 
which we now proceed. 

6. A particular solution for case (B).—Let us suppose that initially (B,); is 
situated to the left of the origin of x and that we are given that 

3 V2 
c2(x) =a? — 


where A,a are positive constants. Hence 


(x<0), (6.1) 


3 1/2 
tar - (ES) 


= — (6.2) 


With these values (5.14) becomes, after evaluation of the integral and simpli- 
fication, 


R(rq) + R( 1g) =2(C— 3Aa (fa? —73)}. 
Without loss of generality we may put C equal to zero and then 


satisfy the equation. Hence, using (5.8), (5.3), the solution of (1.17) corresponding 
to the expression (6.1) for ¢o(x) is 


R(r) = — A(8a-ry, S(s) = | 6.3) 
w= —A(3a—r) 

which implies that 

R_=—A(Za+r,)?, | 6-4) 


From these values, together with (6.2), we find that the second of the equations 
(5.12) becomes 


(— %2)/t, =a, (t2>0), (0.5) 


which gives the position of B, at any instant. Finally (6.4), (5.13) yield 


—4, (6.6) 


which shows that the speed of B, is constant throughout the motion. Moreover, 
by (6.1), this constant is equal to the speed of sound in the far interior of the cloud, 
i.e. in regions remote from the face advancing into the vacuum. 
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It is of interest to find the expressions for u,c in the moving part of the cloud 
in terms of x,t. Using (6.3), equations (5.5) become 


x+(Br+as)t = 
x—(ar+ Bs)t = — 


If we multiply the first equation by (3a+s) and the second by (?a—r) we have 


(6.7) 


($a+s){x+ (Br +as)t} =(2a—r){x —(ar+ Bs)t}, 
which by (1.4), (1.8) becomes after some re-arrangement 
(r+s){x+(a— u)t} =o. 
Again, subtracting the two equations (6.7) we obtain 
(a+ *(3a+s), 
which by (1.3), (1.4), (1.8) becomes 
(r+s)[§t— =o. (6.9) 


We solve these equations by equating the second factors to zero in each case and 
therefore 


(6.10) 


Comparing these with Burgers’ values (4.6) for a homogeneous cloud in which 
k= §, we observe that our values approximate in form to Burgers’ for large values 
of t. 

Turning now to the determination of the speed and position of the advancing 
face B,, we have by (5.10), (6.3) 


— = — A(Za— = —8A(3a—9,)*. 


If we use the first of the equations (6.10) with u, =v, we have also 


A solution of these two equations which, as we shall presently see, corresponds to 
a physically plausible initial state for the expansion is 


% = 3al,— 4 (=) 


1/4 


Finally, let us consider the initial instant of the expansion, which we again define 
by (B,);=(B,); at time ¢;. Instead of equations (4.9), (4.10) we have 


=(— 51) = = = 


(72); =(S2); = — (uz); =0, 


(6.12) 


(6.13) 
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where we have written c; for the value of co(x) calculated at the initial position 
(x,);=(%g) But (6.5) now gives t;=(—.x,);/a=(—-x,),/a, so that the first of 


equations (6.10) is 


Hence also, by (6.11), (v,);=0 and it therefore follows by (6.13) that c;=0, a result 

which we can interpret by saying that initially the face of the cloud is at the absolute - 
zero of temperature.* The second of the equations (6.10), (6.12), respectively, 

when evaluated at the initial instant, are identities, provided that t,=8A/(27a‘), 

and by the first of the equations (6.12) we find 


(*1);= (2); = — 8A 


Thus the constant A in the solution (6.3) depends on the initial position of the 
face of the cloud. 

Summarizing the conclusions which we draw from this example, we may say 
that, if the expansion begins with B, coincident with B,, the face of the cloud 
is then at the absolute zero of temperature and the motion starts with B, moving 
off to the left at the speed a, whilst B, remains momentarily at rest. In common 
with all the solutions of type (5.3), the advancing face is permanently the seat 
of “‘shock”’ conditions. Moreover (6.12) indicates that the high speed of the 
advancing face found by Burgers—three times the velocity of sound in the un- 
disturbed gas of a homogeneous cloud—is unlikely to be realized in nature. 
In our non-homogeneous cloud, with its ‘‘cold”’ initial edge, such a high speed 
would only be attained as a limit after the expansion had proceeded for an infinite 
time. 


Queen Mary College, 
Mile End Road, 
London, E.1 : 
1950 January 29. 


* Incidentally we should then also have at the face a zero root mean square velocity which would’: 
thus be identical with (v,);. 
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Summary 


The paper solves by Riemann’s method the one-dimensional problem of 
the expansion of a gas cloud into a vacuum. The cloud is initially at rest, and 
is homogeneous except for a boundary layer of arbitrary thickness. Motion 
starts in the non-homogeneous layer, and spreads into the cloud. 

The solution uses Riemann’s variables r= 4(3c+ u), s=4(3c—u), where c is 
the local velocity of sound, u the gas velocity. In the first stage of the motion 
there are three moving boundaries. In x>x,(t), there is a vacuum; in 
x_(t)<x<x,(t), r and s are both variable; in x;(t)<x<x,(t), r is constant; in 
x<x,(t), the gas is homogeneous and at rest. The boundaries x= x,(t) and 
x =x,(t) move to the right; and, as x2 increases more rapidly than x,, the two 
boundaries coincide at the end of the first stage of the motion. The boundary 
x=x,(t) moves into the cloud with the local velocity of sound. 

In the second stage of the motion, the quantity 7 is constant everywhere. 
The cloud still advances into the vacuum, and the motion spreads as before into 
the cloud. This state is similar to one discovered by Burgers. The difficulty 
regarding the discontinuity in the initial state of Burgers’ solution is shown to 
be due to the fact that his solution is a limiting case of the present solution 
when the thickness of the initial layer of non-homogeneity is made to vanish. 


1. Introduction.—his paper is concerned with the one-dimensional problem 
of the expansion of a gas cloud into a vacuum. ‘The cloud extends to infinity 
on one side of its infinite plane face; on the other side of the face there isa 
vacuum into which the cloud is advancing. It is assumed that the pressure and 
density satisfy the adiabatic law p =kp’ everywhere; and the effects of viscosity, 
heat conduction and radiation are neglected. 
: The first solution of a problem of this type was obtained by J. M. Burgers.* 
; In this solution, only a layer of the cloud is in motion, the rest being stationary 
; with uniform pressure and density. ‘The boundary between the stationary and 
H moving parts of the cloud moves into the cloud with constant speed a, the velocity i 
; of sound in the gas at rest, whereas the face of the cloud advances into the vacuum 
with speed 2a (y—1). ‘The surprisingly high speed of advance—3a for a cloud 
of monatomic hydrogen with y =3—is a somewhat unsatisfactory feature of this 
solution. A second is that the motion would appear to be due to an unusual type 
¢ of discontinuity at the face of the cloud at the initial instant when the two moving 
‘ interfaces coincide and the cloud is at rest. 
A second solution has been found by G. C. McVittie + for the case y =§, by 
means of Riemann’s theory of the propagation of sound waves of finite amplitude. 


* J. M. Burgers, K. Ned. Akad. v. Wet., 49, 588-607, 1946. 

t Lecture given at the Joint Colloquium on Problems of Motions of Gaseous Masses of Cosmical 
Dimensions held at the Institut d’Astrophysique, Paris, 1949 August 16-19. G. C. McVittie, 
M.N., 110, 224, 1950. 
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McVittie’s cloud is, in some respects, of the same general type as Burgers’; as 
the cloud advances into the vacuum, the depth of the layer of gas in motion increases. 
The speed of advance of the cloud face is, however, no longer constant; Burgers’ 
velocity 3a is the terminal velocity of advance of McVittie’s cloud. 

Despite the elegance of McVittie’s solution, it is not a very satisfactory model. 
He assumes that the stationary part of the cloud is non-homogeneous ; but if the 
adiabatic law held throughout the cloud, this would imply that the whole cloud is 
in motion. This difficulty can be overcome by assuming that the temperature 
of the cloud is proportional to p-! when the gas is at rest but proportional to 

p’-' when the gas is moving*; but this seems a very artificial solution of the 
difficulty. 

A third type of solution is given in the present paper. ‘The cloud is assumed 
to be initially at rest and is homogeneous except in a layer of arbitrary depth near 
the face. This layer starts moving and the face of the cloud advances into the 
vacuum. But in the cloud there are now two moving boundaries, one moving 
to the right, the other to the left; between them Riemann’s variable r remains 
constant. After a certain finite time the former boundary catches up the face of 
the cloud. Thereafter the face of the cloud advances with constant speed 3a, 
and Riemann’s variable r is everywhere constant; in short, the motion is now 
very like that in Burgers’ model. In the second stage of the motion the difficulty 
of the high velocity of the face of the cloud remains; but the time taken to reach 
the second stage can be increased indefinitely by increasing the depth of the 
initial non-homogeneous layer. ‘The Burgers model, on the other hand, corre- 
sponds to the limit as the depth of the initial layer tends to zero, and this explains 
the difficulty of the initial discontinuity at the face of the cloud. 

2. The equations of motion and the boundary conditions.—The motion being 
in one dimension, the dynamical equation is 


in the usual notation. Since, by hypothesis, the adiabatic condition p=kp” is 
true everywhere, we may replace the variable p by the local velocity of sound c 
defined by the equation 


= kyp’-, 


The motion is then governed by the pair of equations 
2c ac 

y-10x 

y—I 0x Ox 


0, (2.1) 


= 0. (2.2) 


* Professor McVittie made this suggestion to me in a letter. I owe him my hearty thanks for 
interesting me in this problem and for his valuable criticisms of my work. 


| 

whilst the equation of continuity is 
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If, in the initial state, the gas is at rest but is non-homogeneous, we have 
u=0, €=¢)(x) when t=0. Equation (2.1) then gives 


so that the initial acceleration is non-zero, and the cloud starts moving. This 
shows that McVittie’s solution is, as stated in the introduction, incompatible with 
the adiabatic law. The equation also shows that, if the cloud lies initially in the 
region x<h say, it is necessary, if expansion occurs, that c;(x) be negative at least 
in some neighbourhood of the face of the cloud. We shall assume that c¢(x) is 
negative in the whole layer of initial non-homogeneity. 

We shall, moreover, consider only continuous initial states of the cloud, so 
that when ¢=0, u=o and p and ¢ are continuous functions of x. It follows that, 
in the subsequent motion, shock conditions cannot occur in the cloud. The face 
of the cloud must, however, be a moving surface of discontinuity. Now if 
x=€(t) is a moving surface of discontinuity, the conditions of conservation of 
mass, momentum and energy require the continuity of 


p(u—£),  pu(u—£)+p and 


across the surface, e being the internal energy per unit mass. It follows in the 
adiabatic case that, since p is zero in the vacuum on one side of the face of the 
cloud, then, in the cloud p and ¢ must tend to the limit zero as we move up to the 
face. No information is obtained concerning the limit of u or the value of €. 

3. Riemann’s form of the equations of motion.—As McVittie has pointed out, it 
is best in this sort of problem to use Riemann’s variables 


r=—+- 


y-I 2 
The equations (2.1) and (2.2) then become 


(3.1) 


(3.2) 
where B=(y—3). 


If one of the variables r,s is a constant, one of these equations is an identity, and 
the other determines the remaining variable. In the general case the equations 
can be linearized by taking r and s as independent variables. We then have 


x—(ar+ Bs)t = ‘ (3-3) 


x+(as+Br)t = — =. (3.4) 


drds as} (3-5) 


From now on we shall consider only the simplest case in which 


y=, = —§. 


c u c u 
3 
or or 

(ar+ Bs) = = 0, 

ot Ox 
. 
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‘The partial differential equation (3.5) simplifies to 
aw 
Ords  (r+s) 


Or > Os = 0, 
of which the solution is 

mee R(r) + S(s) 
r+s 


where R(r) and S(s) are arbitrary functions of r and s respectively. Equations 
(3-1) and (3.2) then give 
R R+S 
§s)t = rts (rts? (3-6) 
Ss’ R+S 


(3-7) 
where accents denote derivatives with respect to the relevant variables. When 
R and S have been determined from the initial conditions, these equations will 
give r and s in terms of x and ¢. 

4. The initial value problem.—Let us consider a gas cloud in which the adiabatic 
law p =kp®* holds everywhere. Initially the cloud is at rest in the region x <h; 
the region x>hA is a vacuum. The portion of the cloud in which x<o is 
homogeneous, whereas the layer 0<x<A is non-homogeneous. The initial 
conditions are therefore 


c=a (x<o), 
(o<x<h), 


where, in view of the continuity assumptions made in Section 2, ¢(0)=a, 
£9(4)=0. In terms of the variables r and s, these conditions are 


r=s=b (x<o), 
r=s=r9(x) (0<x<h), 
where r(0)=b, ro(h)=0, b=§a. 

We assume, in addition, that r,(x)<0, so that each particle of the non-homo- 
geneous layer starts to move to the right. At time ¢ the position of the face of 
the cloud will be x =.x,(t); and on this face r+s=o. 

It is known from Riemann’s work on the propagation of sound waves of finite 
amplitude that a region in which r is constant moves to the right, whereas one in 
which s is constant moves to the left. Hence the initial boundary of the non- 
homogeneous layer will split into two boundaries, and the state at time ¢ will be as 
follows during the first stage of the motion :— 

in the region x<x,(t), r=s=b; 

in the region x,(t)<x<x,(t), r=b, s=s'(x,t); 

in the region x,(t)<x<x,(t), r=r(x, t), $=s(x, t). 
‘The positions of the two boundaries in the cloud are, a priori, unknown; all we 
are given is that r and s are continuous across the boundaries, so that 
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We consider first the region x,(t)<x<x,(t). When t=0, this is the non- 
homogeneous layer in which r=s=r,(x). Hence, by (3.6), (3.7), 


R'(ro) _ + S(70) 


219 475 
S'(r9) R(ro) + S(ro) 
219 


Equating the two values of x, we have 


Ri(ro) + S'(ro) R(to) + Slo) _ 


To 


But as ry is an arbitrary function of the variable x, this implies 
+ =2Ar, 
where A is a constant. We may therefore write 
R(r)=Ar+4(r), S(r) =Ar—¢(r), 
and equations (3.6) and (3.7) then become 


£2. 
r+s (r+s) 
_ , 46s) 


x+(js— §r)t = (4.2) 


(4.1) 


x—(3r— §s)t 


When ¢t =0, either of these equations gives 


$'("o) 


x= 


(4.3) 


where r=r9(x) is the common initial value of r and s. Now, by hypothesis, 
decreases steadily in 0 <x and so has a one-valued inverse function 
x(r9), defined in the range 0<rg<b. It follows that 


=2 rx(r) dr. 


Thus ¢(r) is defined by the initial conditions in the range o<r <6; we have next 
to continue the definition to negative values of r. We do this by considering 
the behaviour of the solution (4.1) and (4.2) as x>~x,(t). 
As x->x,(t), r+s->0. Hence on the face of the cloud, 
r=r,(t)=r(x,,t)>0, $=5,(t) =s(x,, t)<0, 
where 7; +5,=0. Now, if we subtract (4.1) from (4.2) we obtain 
Rt =2 = 
(r+s)° (r+s)? 


(4.4) 


The expression on the right-hand side of this equation must remain finite as 
1-1, S->—1r,. But if we write the equation in the form 


3t = 2{$(r) — 4(s)} —(r +s) $(5)} 
(r+s)* 


si, 
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we see that the right-hand side tends to a finite limit only if ¢(r)=¢(—r). Hence 
¢(r) has to be continued into negative values of r as an even function of r. The 
function ¢(r) is now defined in the range —b <r <b, and we assume that all the 
derivatives which occur are continuous. 

Having defined ¢(r) for negative values of 7, we can find the position of the 
face of the cloud at any instant in the first stage of the motion. Equation (4.1). 
can be written as 


4 _ (r+5)$'(r)- $7) + 
— — §s)t (r+s)? 


and the limiting form of this relation as r—r,, s-> —r, is easily found by 
lHospital’s rule. ‘The result is 


x, —2r,t=44"(r)). 
Similarly, from (4.4) we have 


Hence the position x =x,(t) at the face of the cloud is given by 
(4-5) 
= — 4t. (4-6). 


where 1, satisfies 


From (4.5) we have 
> + 2r,t + => 


But on the face of the cloud c vanishes, and so 2r,=u(x,,t). ‘Thus the cloud 
advances into the vacuum with the speed of the gas at the boundary, a result due to 
McVittie. 

There are several points to note about equation (4.5). In the first stage of 
the motion, 7, increases continuously from 0 to 6. If (4.6) has a solution in this. 
range for positive values of t, it is necessary to assume that $”(r)<o ino<r<b, 
a condition which imposes an additional restriction on the initial data. Again, it 
may happen that (4.6) has more than one positive root ; since x, =2r,, the relevant 
root is the greatest positive root. Other roots are of no importance; for if a 
value of r inside the cloud satisfied the relation ¢’(r) = — 4t, it would not follow 
that the corresponding value of s was —r. 

The inner boundary of the outer layer is fixed by the conditionr=b. Hence, 
if we write x =x,(t), r=b, s =s,(t) in (4.2) and (4.4), we find that this boundary is. 


given by 
where * 
In the region x,(t)<x<x,(t), r=6; and s is given by 


* We note that ss=6 implies t=o, thereby verifying our assertion that the inner boundary of the- 
initial non-homogeneous layer splits into two. 


| 
| 
| 
i 
| 
| 


-constant everywhere. 


-(4.5) and (4.7) equal values of x,(7) and x,(7) when r,=b, s,=—6. This is, in 


at the instant 7, and must therefore be the same function as in (4.9). Hence, 
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which has the solution 
x + (3s — §b)t =g(s). 


The function g(s) has to be chosen so that s is continuous across x =x. Comparing 
this equation with (4.7), we see that 


The third moving boundary x =.,(t) is specified by s=5. This gives 
2b 


But, by (4.3), ¢(6)=0, and so x,;=— This boundary has velocity 
x, = — 3b= —a, where a is the velocity of sound in the gas at rest. 

5. The second stage of the motion.—The motion described in Section 4 will 
persist so long as the three boundaries are distinct. But if the two boundaries 
x=x,(t) and x= x,(t) coincided at an instant 7, at that instant we should have a 
moving layer x,(7')<*<.x,(T) in whichr is constant; and we could take the state 
at the instant 7 as a new initial state. Inthe subsequent motion, r would remain 


Since r+s=0 on x=x,(t) and r=6 on x=x,(t), the two boundaries would 
coincide at the instant t= T only if (4.6) and (4.8) gave the same value of T and 


fact, the case; the values are 


T=- 
x(T) =x,(T) =2bT + 14"(b). 
Moreover, at this instant, the two coinciding boundaries have the same velocity 


2b =3a. 
When t>T, r will be constant everywhere, and s will be given by an equation 


x + (js — =g(s) 
inthe moving partofthe gas. ‘The function g(s) has to be chosen to fit the conditions 


when t> 7, the motion is confined to a layer x,(t)<x<.x,(t) and is specified by 
r=b, 


x+(3s—§b)t = (as? 
The boundary x =.,(t), corresponding to s=4, is given by x3 = — as before. 
The face of the cloud x = x,(¢) is no longer given by (4.5), but is obtained by making 
.s—> —b in (5.1); the result is 


=2bt + (5.2) 


Since 


the face of the cloud advances with a speed equal to three times the velocity of 
sound in the gas at rest. The model suffers, then, from the same defect as 


x 
43 
| 
| 
ott 
| 
| 
i q 
x, =2b = 3a, 
Pang 
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the depth of the initial layer of non-homogeneity. 


V2 
=b (: i) 
It follows from (4.3) that 


layer is given by 


At the face of the cloud x =.x,(t), we have r= —s=r,(t). 
r, =6*t 3h, and hence, by (6.1), 
xy =h + Sh 
The velocity of the face of the cloud is «, = 267t/3h in the first stage of the motion, 
which lasts, as we shall see, until x, attains the value 26. 


At the inner boundary of the outer layer x=x,(t), we have r=b, s=s,(t). 
Hence 


2b*t 
s,(t)=b — —- 
) 3h 
and 
2bt 2b?t? 
=—+—;. 6. 
The velocity of the inner boundary is 
2b 
=— +=, 
so that the inner boundary starts to move with the local velocity of sound 26/3. 


Since 


the inner and outer boundaries of the outer layer coincide at the instant t = 3h/b. 
At this instant, x, =x,=4h, x, =x, =2b. 
In the region x,(t)<x<x,(t), we have r=b, and hence 


x + (4s— §b)t =g(s), 
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Burgers’ model in respect of this high speed of expansion, but with the important 
difference that the time T taken to attain this speed can be increased by increasing 


6. A simple example.—A very simple example is provided by the case when. 


when 0 <r <6, no constant of integration being necessary, since it would disappear 
from (4.1) and (4.2) if it were inserted. As ¢(r) is already an even function, 
the definition holdsin —b<r<b. It readily follows that the motion in the surface 


(6.2) 


From (6.2), we have 


hr* 
=hr*? — 
r-s=—. 
| 
| 
| 
| 
| } 
= bt \2 
18 
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where g(s) has to be chosen so that s=s, when x=. It is easy to show that 


x+(4s—$b)t = (6.4) 


But on x=x,(t), we have s=6, and hence 
= — 


just as in the general case. ‘This completes the description of the first stage of 
the motion. It lasts until the instant t = 34/6; hence the time taken for the outer 
boundary to attain the velocity 3a can be made as long as we please by making h 
large enough. 


In the second stage of the motion, r=5; and s is given by 


x + (4s— 3b)t = 


On the face of the cloud we now have s = —), and so 
x, =2bt — 2h, 
which advances with constant speed 2b, whereas on the inner boundary 
= — 
In conclusion we observe that, if we make h—>0, we find from (6.4) that 
x + (4s— =o. 


This is Burgers’ solution. We have derived it by considering an initial state, in 
which the cloud is at rest in x <h, defined by 


r=s=b(x-3) (0o<x<h), 


r=s=b (x <0), 


and then making A tend to zero. This explains the singularity in the initial state 
of Burgers’ solution to which McVittie refers. 
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THE CONTINUOUS RADIATIVE ABSORPTION 
CROSS-SECTION OF SINGLY IONIZED POTASSIUM 


M. F. Seaton 
(Communicated by H. S. W. Massey) 
(Received 1950 March 23) 


Summary 

After obtaining exchange wave functions of K III the continuous 
absorption cross-section of K II is calculated over the energy range o to 
13°53 e.V. of the ejected electron, using both the dipole length and the dipole 
velocity formulae and allowing for the effect of electron exchange. The 
recommended mean value is taken to be x 10~'"v,/y cm.*®, where 
is the frequency at the spectral head. 

The cross-section data are used to evaluate the “‘ modified dilution 
factor” introduced by Strémgren in considering the ionization balance in 
interstellar space. The possibility of captures to excited states introduces a 
correction factor of 430 (T,=10*), 290 (T= 10°), 160 (T,=10*), where T, is 
the electron temperature (deg. K.). The abundance of interstellar potassium 
is briefly discussed. 


I. Introduction 
While the quantal theory of continuous radiative absorption cross-sections 
has been applied to a fair number of neutral atoms, comparatively little work 
has been done on positive ions.* Singly ionized potassium is a case of some 
interest as representative of a group of positive ions with outer p* shells which 
should be amenable to comparatively accurate calculations. It is also of 
importance in the study of the ionization balance in interstellar space. The 
continuous absorption cross-section of KI defies accurate calculation due to great 
sensitivity to details of the wave functions, and it is therefore fortunate that 
experimental results are available for this case, enabling a calculation of potassium 
abundances to be made for HI regions.t For the extension of these calculations 
to HII regions it is necessary to have data on the absorption cross-section of K II. 


Il. The Absorption Cross-section 
(1) Quantal Formulae for the Cross-section.—After substituting for the 
fundamental constants involved, the quantal formula for the continuous 
absorption cross-section {| may be put in the form 


a, =8:55X 10 +€)@, {C1 + of, 1} (1) 
where the ionization potential J and the energy of ejection « are both expressed 


in Hartree atomic units (13:53 e.V.). @,, a factor allowing for the distortion of 
the wave functions of the passive electrons, is given by 


nl (all passive electrons) 


0 


* For a general review see D. R. Bates, M.N., 106, 432, 1946. References to later work are 
given by D. R. Bates and M. J. Seaton, M.N., 109, 698, 1949. 
+ For the experimental results for KI see R. W. Ditchburn, J. Tunstead and J. G. Yates, 


Proc. Roy. Soc. A, 181, 386, 1943, and for a theoretical discussion, D. R. Bates, Proc. Roy. Soc. A, 
188, 350, 1946. 


1D. R. Bates (1946). 
18* 


& 


| 
i | 
| 
it 
| 


248 M. #. Seaton, The continuous radiative Vol. 110 


where P*(nl|r) rand P{(nl|r)/r are the radial wave functions, normalized to unity, 
of the passive electrons in the initial and final states. The factors C,_, and C,,,, 
which have been tabulated by Bates (1946), arise in averaging over the initial 
state, summing over the final state, and integrating over the angular coordinates. 
The fundamental matrix elements, o,,,, can be expressed in terms of the dipole 
length (o,,,(L)) or in terms of the dipole velocity * (,,,(V)). In the former 
case we have 


0141(L) = P\(nl|r)P{e, 1+ 1\|r)r dr (3) 


and in the latter 


where P'(nl|r)/r and P(,l+1|r)/r are the radial wave functions of the active 
electron in the initial and final states, P'(nl|r) being normalized to unity and 
P%e,1+1|r) to asymptotic amplitude 1/e'*. As remarked by Bates and Seaton 
(1949) the two expressions for the matrix elements can only be expected to give 
identical results if the wave functions are solutions of the Schrédinger equation, 
while the wave functions with which we are concerned are the best separable 
wave functions for the determination of the energy. 

(2) Hartree-Fock Wave Functions of K 111.—The equations satisfied by the 
free electron wave functions of K II involve the wave functions of the K III core. 
Although wave functions were available with and without exchange for K II, 
none appears to be available for KIII. However, wave functions (without 
exchange) have been calculated for the iso-electronic case of Cl and Cl.{ 
Hartree, Kronig and Petersen have remarked that for these the inner part of the 
charge distribution remains almost unaffected in going from Cl to Cl (due to 
the simultaneous diminution and contraction of the total charge cloud). Thus 
the tabulated functions P(1s), P(2s) and P(2p) are identical for Cl and Cl’. An 
estimate of the wave functions of KIII was therefore made in the following 
manner :— 

(i) P(1s), P(2s) and P(2p) of K III were taken to be the same as for KII. 

(ii) For P(3s) the function p= P(3s|Cl)— — P(3s|CT), which is nowhere 
greater than 0-009, was plotted against p =r/r,,, where 7,, is the principal maximum 
in the chlorine functions. The value of p on the p scale of P(3s| K II) was then 
read off and the function P(3s|KII1)=P(3s|K11)+p was calculated and 
normalized as the required estimate. 

(iii) It was not possible to use the same procedure for the P(3p) functions, 
since P(3p|Cl-) is much too diffuse with exchange neglected (cf. Hartree, 1938). 
The procedure adopted was to calculate p= P(3p| Cl) — P(3p| Cl-) as a function 


of p and then to put 
P(3p| KU) 
P(3p| Cl) 


*S. Chandrasekhar, Ap. 7., 100, 176, 1944; 102, 223 and 395, 1945. S. Chandrasekhar and 
F. H. Breen, Ap. 7., 104, 430, 1946. 

+ D. R. and W. Hartree, Proc. Roy. Soc. A, 166, 450, 1938; and D. R. Hartree, Proc. Roy. Soc. A, 
143, 506, 1934. 

} D. R. Hartree, Proc. Roy. Soc. A, 141, 282, 1933; and D. R. Hartree, R. de L. Kronig and. 
H. Petersen, Physica, 1, 895, 1934. 
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so that p’ tends to zero as P(3p| K II) tends to zero. The function 
P(3p|K III) = P(3p| K II) +p’ 


was calculated and normalized as the required estimate. 

The above estimates of the Is, 2s, 2p and 3s functions were used without 
further improvement. The estimate of P(3p|KIII) was used as a first 
approximation in the solution of the Hartree-Fock equations for KIII. The 
final wave function obtained, which is given in Table I, was found to be quite 
close to the original estimate (nowhere differing from it by more than 0-007). 


TABLE I 
P(3p) for K III, with exchange 
€3p3p= 3°43 A.U. 


“333 


"252 
“000 
132 
"259 


ry is in the usual atomic units, 0°528,; x 10-8 cm. 


As a check on the wave functions the ionization potential of KI1I was 
calculated. - Neglecting the variation in all wave functions except P(3p), this is 
given in atomic units by 


I=[—20F(3p3p) + 2 — 5€spap) 


—[—30F ((3p3p) + F2(3P3p) — %€spap (5) 
where 


Finint) = { P%nl\r’) (Fy ars P%(nl|r’) dr, (6) 


r’ 


& 
ie 


r r=0 
0-0 0°000 —o'582 
‘006 — 735 
+023 — +839 
03 ‘047 — 
"04 ‘076 I'l 
"05 108 1°2 — ‘919 
-06 142 
—o'856 
0:08 0°210; 1°6 j 
"10 2745 1°8 — 640 
“12 2-0 — +528 | 
14 "375 2°2 428 
18 "429 2°6 — +268 
*20 2°8 — +209 
"22 “439 — +161 
"24 "429 — 
‘26 “411 — *094 
+28 +385 3°6 — ‘O71 
4°0 — +040 
0°35 ° — 
40 + 4°5 —o-o18, 
45 50 — 968 
"50 5°5 — "004 | 
55 6-0 — "002 
“60 6°5 — ‘oor 
| 
| 
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and the e’s are the energy parameters in the Fock equations. The expression 
for J was found to be very sensitive to the details of the wave function P(3p). 
The final value obtained, /=2-47A.U. was estimated to be correct to within 
+5 per cent so far as the self-consistent procedure is concerned. ‘This is to be 
compared with the experimental value J=2-341a.u.* The experimental value 
of J was used in all formulae. 

The factor @,, was found to be 0-979. 

(3) Matrix Elements and Cross-section at the Spectral Head.—The free electron 
wave functions were calculated at the spectral head from the Hartree-Fock 
equations using the usual method of successive approximations in which the 
first approximation is obtained by neglecting the exchange terms. The 
normalization was carried out by a method due to Strémgren (cf. Bates and 
Seaton, 1949). The matrix elements, calculated using all combinations of 
exchange and non-exchange wave functions, are given in Table II. We denote 
a wave function with exchange by P,, and without exchange by Py. The values 
shown in ordinary type are obtained from the dipole length formula and those 
in italic type from the dipole velocity formula. The spectral head cross-section 
is given in terms of the o’s by 


@, = 0°392(a7_, +207, ,)10-!7 cm.?. 


Taste II 


of, , for €=0 


2 


P (es) P,(ed) P (ed) 


P(3p) 0°67, 0°34 || Po(3p) 3°8, 2:0 5°9, 2°3 


P(3p) 0°52, 0°38 || PX3p) | 2°5 2°5° 


It is apparent from Table II that the contribution from p—s transitions is 
much smaller than that from p—-d transitions. This is a common characteristic 
of systems with outer p* shells. The function P(es) is small and oscillatory 
in the region giving the major contribution, while P(ed) has its first node 
comparatively far out and rises to a large maximum in the region where P(3p) is 
large. The exact value of o7_, is of little importance in calculating the cross-section 
and will not be considered further. 

For the d wave it is instructive to compare the position with that of the neutral 
atoms with outer 2p electrons, of which oxygen is a typical example t : 

(i) Sensitivity to details of the atomic field.—For oxygen the radial equation 
for the d wave is so dominated by the term — 6/r? in the region where the potential 
differs much from its asymptotic form that details of the atomic field are not 
important. In consequence, there is very little difference between P,(ed) and 
P (ed) and a hydrogenic wave function gives a very good approximation. Due 


* J. C. Boyce, Rev. Mod. Phys., 13, 1, 1941. 
+ D. R. Bates, M.N., 100, 25, 1939. Also D. R. Bates and M. J. Seaton (1949). 
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to the greater nuclear charge the term —6/r? is not so important for KII and 
the d wave becomes quite sensitive to details of the atomic field. Exchange 
appreciably modifies P(ed), the effect being repulsive, and the hydrogenic 
approximation is rather poor. 

(ii) Cancellation effects.—For oxygen the initial node of P(ed) is so far out 
(r =9°7A.U.) that no significant cancellation occurs in the integral for the matrix 
element. The initial node of P(ed) for KII occurs at much smaller radial 
distances (r = 3-6 without exchange, r = 4-5 with exchange) but it is still sufficiently 
far out to exclude any serious cancellation effects.* The contraction is mainly 
due to the fact that we are now dealing with a positive ion and only to a lesser 
extent due to the increased nuclear charge. 

(iii) Sensitivity to details of the bound electron wave function.—The 
sensitivity of 07, , to exchange in P(3p) is surprisingly small. This is due to the 
comparative compactness of the integrand, a significant contribution coming 
from both sides of the radial distance at which P,(3p) and P,(3p) cross over. 


INTEGRAND 


F, (3p 3p) 
INTEGRAND 


A 


2 3 4 
Fic. 1.—Integrands of dipole length, dipole velocity and energy integrals. 


There is a marked discrepancy between the dipole length and the dipole 
velocity formulae 7, and it is of some importance to decide which is likely to be 
the more accurate. In Fig. 1 we have plotted the integrands of o,,,(L) and 


* A convenient measure of the degree of cancellation in the integral is given by the quantity 
(1—D), where D is the ratio of the smaller to the larger of the contributions of opposite sign. For 
p — d transitions (1 — D) is in all cases greater than 0-8. For comparison we note that Bates found 
(1—D) for K I (s > p transition) to range from 0-388 to 0-015 according to the approximation used. 
For oxygen, on the other hand, (1—D) is never appreciably below unity for p > d transitions. 

+ That this is a good deal greater than in the case of oxygen is not unexpected in view of the 
greater complexity of the system. 


£+1(V) 
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0;4,(V) and also the integrand of a typical energy integral, F)(3p3p) 
(cf. equation (6)). It will be noted that the dipole velocity integrand is large in 
the same regions as the energy integrand, whereas that of the dipole length is 
still large for radial distances where there is a negligible contribution to the 
energy. Bearing in mind the principle on which the wave functions are 
calculated, this would certainly suggest that greater reliance is to be given to the 
dipole velocity integral, and the smaller sensitivity to details of the wave functions 
would appear to substantiate this view. It does not, however, seem possible to 
reach any final decision on this point until a comparison of theory and experiment 
can be made for at least certain representative substances.* In the meantime 
the wisest procedure is probably to take a simple mean in order to reduce the 
maximum possible error. This gives a, = 3:2 < 10-'? cm.? at the spectral head for 
the case of full inclusion of exchange. 

(4) Frequency Variation of the Cross-section.—In order to calculate the 
frequency variation of a, the function P,(ed) was calculated for «=0-0, 0-2, 
0-4, 0°6, 0-8 and 1-0a.u. The numerical work was considerably reduced by 
neglecting the frequency variation of the exchange terms, a procedure which 
is justified by the fact that these terms are only large in the region where the 
frequency variation of P,(ed) is not excessive and by the lack of sensitivity of 
P,{ed) to details of the exchange terms. 

The function P,(es) was calculated for «=o and for e=1-0. ‘The corre- 
sponding values of o7_,(L) were found to be 0-52 and 0-46 and those of o7_,(V), 
0-38 and 0-13. Intermediate values of o7_, were calculated using a simple linear 
interpolation. 

Fig. 2 gives the frequency variation of a, using the two basic formulae. ‘The 
dipole length cross-section is found to be closely proportional to 1 v while that 
of the dipole velocity falls off a little more rapidly. Consistent with the mean 
value at the spectral head we will take the cross-section to be given as a function 
of frequency by 
Vy 2 
cm.”; (7) 
where vy, is the frequency at the spectral head. 

It is not easy to estimate the reliability of this result. On the one hand, the 
absence of serious cancellation in the matrix elements makes any gross error 
unlikely. On the other hand, the difference between the dipole length and velocity 
formulae would make one hesitate to assert that the final recommended value is 
correct to within less than a factor of two. It will be seen that the uncertainty 


involved is of little importance for the application to the problems of interstellar 
matter. 


@,=3°2 X10°"" 


Ill. Astrophysical Application 
(1) Interstellar Ionization Equilibrium.—'Vhe analysis of interstellar absorption 
lines can usually only be expected to give the abundance of an element in a 
particular stage of ionization; in order to calculate chemical abundances it is 


* There is some indirect evidence available on this point. D. R. Bates and A. Damgaard 
(Phil. Trans. Roy. Soc. A, 242, 101, 1949) consider that it seems likely that the velocity formula gives 
a better value for the oscillator strengths of spectral lines for complex systems. However, for 
bound-bound transitions the integrals are necessarily more compact and the problem would not 
appear to be so acute. (It is of interest to note that Bates and Damgaard find that exchange gives 


a considerable improvement.) Neon would appear to offer the best possibilities for a direct com- 
parison of theory and experiment. 
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a, in units 
10-!7 cm.? 


5 


DIPOLE LENGTH 


DIPOLE VELOCITY 


L 
0-2 0-4 0-6 0-8 1-0 

€ in atomic units (13-53 e.V.) 
Fic. 2.—Frequency variation of the cross-section using the dipole length and 
dipole velocity formulae. 


necessary to study the ionization equilibrium which is set up. In doing this we 


must bear in mind that the conditions are about as far from thermodynamic 
equilibrium as it is possible to imagine *, and in particular we cannot expect the 
principle of detailed balancing to be applicable. Strémgren + has shown that the 
correct procedure is to equate the rate of photo-ionization from the ground state 
of an atom or ion to the sum of the recombination rates to all excited states. This 
procedure gives the ionization equation 


where &, is the electron density, #{, # the densities of atoms or ions in the 
ground state in two consecutive stages of ionization, g},qj the corresponding 
statistical weights. 7; is the radiation temperature, T, the electron temperature 
and / the first ionization potential. The factor W is given by 


(0) 


é. 
where a‘ is the absorption cross-section for the kth excited state of the atom or 
ion concerned (weight function q,), k= 1 corresponding to the ground state. E, is 
the kinetic energy of the electrons and G, is a factor defined by 
2hv® 


(x0) 


* Cf. A. S. Eddington, Internal Constitution of the Stars, Cambridge, 1930. 
+ Bengt Strémgren, Ap. 7., 108, 242, 1948. 
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where J, is the radiation intensity. (When J, is given, the choice of 7; is quite 
arbitrary; in fact 7; could be eliminated from equation (8) altogether.) 

The familiar Eddington-Saha equation has the same form as (8) with W 
replaced by a dilution factor G. For this reason W is referred to by Strémgren 
as the ‘‘ modified dilution factor’. In order to calculate the effect of recombina- 
tions to excited states it is convenient to adopt the following procedure. We first 
coleulete W by ge: all recombinations to excited states, and assuming 
that a! varies as 1/v?. Using the intensity data of Dunham *, as recommended 
by Strémgren, this gives for the ionization of K II 


(kT, 
W(1)= 2: 87 x 


(11) 


Secondly we may calculate W, using the same intensity data, the calculated 
expression for the cross-sections, and the inclusion of the sum over excited states. 
Calling this quantity W(2) the ratio W(1)/W(2) is termed the “ capture 

correction’ by Strémgren. It has been calculated for Nal, K1, Cal and Call I 
and varies from 2:6 to 1600. In calculating the sum over excited states a 
hydrogenic approximation for the absorption cross-sections, as given by Kramers’ 
formula for HeII, was judged to be sufficiently accurate.t ‘The denominator 
of W(2) may then be put in the form 


where the He II sum is taken over the 3d electron and all m equal to or greater than 4. 
After the first few terms it is legitimate to replace the He II sum by an integral 


: which may be evaluated analytically. Putting x, = a 
head frequency of the mth excited state, we have for the higher excited states 
of He Il 


where v, is the spectral 


| evn ( | dx, tere | + +7}, 
0 Ty’ 


y being Euler’s constant (05772). ‘The He II sum for n > 4 is given in Table III. 
It is found that the higher excited states, for which the cross-sections of K II and 
Hell would not be far different, give a very considerable contribution to the 
sum, which justifies using this approximation. 


* Theodore Dunham, Jr., Proc. Amer. Phil. Soc., 81, 277, 1939. The ionization potential of 
; KI being 31°67 e.V. it will be realized thata very considerable extrapolation of black-body intensities 
has been made in using Dunham’s data. Further work on the ultra-violet spectrum of hot stars 
may call for a revision in the value of the quantity W(1) (equation (11)), but the ratio W(1)/W(2) is 
not likely to be greatly affected. 
t See D. R. Bates, R. A. Buckingham, H. S. W. Massey and J. J. Unwin, Proc. Roy. Soc. A, 
170, 322, 1939, for a discussion of this approximation. For the Kramers formula, cf. D. R. Bates 
(1946). The correction factor g was taken to be unity. 
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Tasce III 
n=4 0 
for He II in units a,v{ of He II 
10? 2°32 
108 4°54 
6°8, 


In Table IV we give the “ capture correction’’ for KI1, W(1)/W(2), for 
electron temperature 10, 10° and 10 deg. K., together with the corrected ratio. 
Rew Rev 


TABLE IV 
T, W(1)/W(2) Rg 
107 430 
10° 290 0066 
104 160 0°38 


(2) Potassium Abundances.—Since for HI regions we can assume a complete 
cut-off of radiation beyond the ionization potential of hydrogen the ionization 
of KII is only of importance in HII regions, where we expect the electron 
density to be high (@,=@y). From the values of the ratio iy 
given in Table IV we can see that so long as &, is greater than I per cm.* the 
ionization of K II will not be important. 

Strémgren (1948) has made a very full analysis of abundances as derived 
from the absorption lines in x? Orionis, treating separately the hypotheses that 
the absorbing cloud is an HI region and that it is an HII region. Using the 
present ionization data we may complete his calculations for potassium. The 
results are given in Table V. 


TABLE V 
H II region 
10 7X 1077 0°06 
I 9 x 10-8 0°08 
orl 30 x 10-5 0°30 


It is interesting to note that for all assumed conditions (except an H II region 
of low density, which appears to be improbable) the abundance ratio #,/By,g 
is close to 0:06 for both HI and HII regions.* This is due to the similar 
mechanism involved if ionization of K II be neglected, since ionization of Na II 
(ionization potential = 47-2 volts) is certainly unlikely to occur. 
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VENUS AND THE EARTH’S INNER CORE 


K. E. Bullen 


(Received 1950 February 22) 


Summary 


On the hypothesis that the density changes inside the Earth, down to a 
depth of 5000 km., are essentially due to pressure and contorm to the writer’s 
compressibility-pressure Earth model, it is shown that, if Venus and the 
Earth are of the same primitive composition, it is more probable than not 
that Venus contains an inner core. This supports the evidence previously 
found (on the same hypothesis) from the case of Mars that the Earth’s inner 
core is chemically distinct from the rest of the Earth. A formal relation 
between the rotation period and the ellipticity of Venus is derived and i 
discussed. 


The use of the Earth model * (Model B) based on a compressibility-pressure 
hypothesis leads to the suggestion ¢ that Mars may have an inner core, and it 
therefore becomes desirable to re-examine the constitution of Venus on the 
basis of Model B. ‘The notation in the present paper will be that of Papers I 
and III {; M and R will here denote the mass and radius of Venus. 

2. If Venus has an inner core in the same proportion by mass as in the Earth, 
the mass of this core, determined in the same way as for Mars (Paper III, 
Section 4), is 1:26 x 10?7gm. The mean density of an inner core in Venus can be 
estimated as in the case of Mars, using the equation (38) of Paper II. As a first 
approximation, the mean pressure in this core was taken as 2-5 x 10!" dynes/cm.?, 
leading to an estimated mean density of 14-7 gm.;cm.*. The corresponding first : 
approximation to the radius of the inner core of Venus was 1270 km. 
Applying the density- -pressure figures in Table III of Paper II to Venus, and 
i taking the crustal layers as in Paper I to be equivalent to a single layer of 
thickness 30 km. and density 2-8 gm./cm.*, it is found that on the data 


(i) M=4-88 x 1077 gm.; R=6200km., 
(1) M=4-88 x 1077 gm.; R=6300km., 


the mass inside a sphere of radius 1270km. in Venus is 0-178 x 107, 
0-100 x 10”? gm., respectively; the corresponding mean densities are 20-4, 
11-8 gm./cm.°. 

Interpolating linearly, we find that the data 


(iii) M-=4-88 x 102? gm.; R=6267 km. 


are consistent with the presence of an inner core in Venus of mean density 
14:7 gm./cm.* and radius 1270km. Also, the data 


(iv) M=4-88 x 1077 gm.; R=6312 km. 


* M.N., Geophys. Suppl., 6, 50, 1950; this paper will be subsequently referred to as Paper II. 
t+ M.N., 109, 688, 1950; this paper will be subsequently referred to as Paper IIT. 
t M.N., 109, 457, 1949. 
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are found to be consistent with the absence of an inner core in Venus. Thus. 
(iii) and (iv) are consistent with an identical primitive composition for Venus. 
and the Earth, assuming Model B to be relevant to the Earth and that the density 
change at the outer boundary of the central core is due solely to pressure; 
(iii) assumes that the Earth’s inner core is composed of chemically distinct 
material, while (iv) assumes that the Earth’s inner core is due to the prevailing 
pressure. 

The figures (i) are the most probable values of M and R given by the 
astronomical observations; the estimated standard error of M is 0-03 x 1077 gm., 
and of R, 50km. Thus the data (iii) are practically within the standard errors 
of the observations, while (iv) are somewhat farther away. ‘Thus, as in the case 
of Mars, if the Earth Model B is relevant, it is more probable than not that Venus. 
has an inner core. As shown in Paper I, the use of the Earth Model A*, without 
assuming an inner core, is consistent with the data 


M =4°88 x 1077 gm.; R=6285km., 


which are not seriously different from (iii). But the case of Mars makes a sharp. 
discrimination favouring Model B against Model A, if the terrestrial planets. 
are of uniform primitive composition. 

Thus the new calculation for Venus gives some small additional support to 
the view that the terrestrial planets contain inner cores chemically distinct from. 
the remaining materials. 

The suggestion made in Paper I on the desirability of improved precision in 
observations of the diameter of Venus is still very relevant. If the uncertainty 
could be narrowed to the order of that for Mars, it would be possible, on the 
theories we have been following, to discriminate much more sharply between 
various competing Earth models. 

3. In Tables I and II the density and pressure distributions in Venus are 
shown on the hypotheses that the Earth Model B is relevant, and that Venus has. 
not and has an inner core, respectively. In preparing Table II, a second 
approximation has been made from the figures in the first paragraph of Section 2, 
resulting in a reduction of 0-6gm./cm.* in the mean density of the Venus inner 
core and an increase in the mean radius to 1330km. The value 6267km. 
for R is not sensibly affected in the second approximation. It is likely that if 
Venus has an inner core, there would be a transition layer corresponding to the 
layer F in the Earth, and this would have a slight effect on the density distribution 
in Venus. In view of the limited knowledge concerning the region F, it seems. 
best to ignore this region for the present in constructing Table II; also the 
extent of density increase inside the inner core is fairly uncertain. On the other 
hand, Table I, which is not complicated by an inner core, is not subject to these: 
particular uncertainties. 

4. The values of J) MR? yielded for Venus using (i) the Earth Model A, 
assuming that Venus has no inner core, (ii) the Earth Model A, assuming that 
Venus has an inner core, (iii) the density distribution of Table I, (iv) the density 
distribution of Table II, are 0-360, 0-354, 0-368, 0-362, respectively. 

If e is the ellipticity of Venus, and m the ratio of the centrifugal force at the 
equator to the gravitational attraction, the corresponding values of e/m given on 
the Radau-Darwin hydrostatic theory are 0-172, 0-168, 0-177, 0-173, respectively. 


*K. E. Bullen, Theory of Seismology, Cambridge, 1947, p. 218. 
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Pp 

0007 10!” dynes/cm.?* 
0°025 
0'025 
0'280 
0°54 
0°67 
1°04 
1°25 
1°33 
1°33 
1°43 
1°63 
1°78 
1°97 


p 
0-007 dynes/cm.? 
0°025 
0°025 
0°123 
0°255 
0°39 
0°52 
0°66 
0-92 
1'05 
1°33 
1°33 
10°13 1°58 
1467 10°45 1°84 
1330 10°55 I‘g! 
1330 13°3 1-91 
° 15°3 2°5 


Let the period of rotation of Venus be 20y days. Then 


(1) 


where Q=3-64x10-*y"! radians/sec. We shall take R=6200+100km., 
M =(4-88 + 0-03) x 1077 gm., e/m=0-172+0-005. Substituting in (1), we have a 
formula connecting e with y, namely 


e =(1-68 + 0:07) x 10-8y-?, (2) 


The formula (2) could of course be significant only if y were much less than 
unity, which is not probable on present indications of the rotation period of Venus. 
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I 

d r p 
30 km. 6282 km. 3°32 gm./cm.* 
97 6215 3°36 

97 6215 3°87 
aes 400 5912 4°06 
800 5512 4°25 
1200 5112 4°43 
1600 4712 4°59 
2000 4312 4°74 
2400 3912 4°87 
Bt 2800 3512 4°99 
3200 3112 5°11 
3600 2712 5°21 
4000 2312 5°32 

4236 2076 5°57 

4236 2076 9°74 
4400 1912 9°90 
4800 1512 10°18 
ae 5200 1112 10°38 i 
5600 712 10°54 
6312 ° 10°63 
II 

95 
| 
1600 
2000 
2400 
2800 
3200 
3600 
4020 
4020 
4400 
4800 
4937 
4937 
e RQ 
we 
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If y is of the order of unity, the value of e would not be resolvable by direct 
observation, and would in any case almost certainly be dwarfed by surface 
inequalities sustained by the strength of Venus. Thus on present indications 
observations of the figure of Venus are not likely to be helpful. 

5. The presence of a chemically distinct inner core in the Earth, Mars and 
Venus is now seen to be independently supported by observations of 


(a) the mass and radius of Mars, 
(6) the figure of Mars, 
(c) the mass and radius of Venus, 


assuming that Model B fits the Earth. The cautions attaching to the use of 
Model B have already been stated in Papers II and III, and the same need 
remains for further tests that will assist in discriminating between Model B and, 


for example, Ramsey’s proposed modification * of Model A, before too definite 
conclusions can be drawn. 


University of Sydney : 
1950 February 14. 


* M.N., 108, 406, 1948; M.N., Geophys. Suppl., 5, 409, 1949. 
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OCCULTATIONS OF STARS BY THE MOON OBSERVED 
THE NIZAMIAH OBSERVATORY, HYDERABAD, 
DURING THE YEAR 1949 


(Communicated by the Director) 


(Received 1950 March 8) 


TABLE I 


yA 


I 
2 
3 
4 
5 
6 
7 
8 


” 


The last column indicates the following observers : 
M.G. : Mahomed Ghouse. 
C.H.: C. Hanumantha Rao. 


The instrument used throughout was the 15-inch Grubb refractor. The- 
corrections to the standard clock were determined from observations with the- 
transit instrument. The reductions were carried out by S. Aravamuthan 
assisted by C. Hanumantha Rao. No correction was applied to the observed 
time. 


: || Star Mag. Phase Obs. 

B.D.—5 6048 7°4 jon. 16 7 64 M.G. 
B.D.—5° 6056 "7 5 16 41 27-0 M.G. * 
88 Pisc 6:2 7 17 44 42°0 M.G. 
B.D.+6° 185 8-3 7 18 6 20°5 M.G. 
B.D.+ 10° 257 7°4 8 13 24 24°3 M.G. 4 
oS 268 Arie 61 a 8 17 49 57°3 M.G. 
107 B Auri 6°5 12 19 1 M.G. 
49 Auri 13 I5 51 13°0 M.G. 
9 Libr 4°7 a 22 22 44 41°8 M.G. { 
10 B.D.—7° 6036 6-4 Feb. 1 14 23 17°7 M.G. 
II 8 Arie 4°5 6 13 41 43°2 
12 36 Taur 5°7 7 14 36 2971 C.H. 
13 B.D. + 26° 775 6-9 8 14 16 52°2 C.H. 
14 B.D.+27° 1337 10 18 23 M.G. 
15 B.D. + 25° 731 7°5 Mar. 7 17 41 48°1 M.G. 
16 25 Gemi 6°5 9 15 28 58-9 M.G. 
17 B.D.—26° 1656 6:8 10 18 56 57°3 M.G. 
18 B.D.+24° 655 9°0 April 3 14 18 53°6 
19 B.D. 1920 71 7 18 45°4 C.H. 
20 28 Canc 6°1 7 18 56 25°5 C.H. 
21 B.D.+20° 2318 6-7 8 18 25 50°5 C.H. 
22 B.D. + 10° 2230 7:0 May 7 19 29 45°4 C.H. 
23 B.D.—13° 3824 6°7 June 7 91 24°97 C.H. 
‘S 24 44 Pisc 6:0 Nov. 2 18 32 30°6 M.G. . 
| 25 B.D.—25° 14589 72 24 14 58 56°7 M.G. 
a 26 B.D.— 12° 6327 6:8 re 27 13 25 30°0 M.G. 
i 27 B.D.—11° 5912 6:7 27 17 35 2°1 M.G. 

a 28 50 Agar 5°9 Dec. 24 14 39 54°4 M.G. t 
29 23 537 68 31 20 50 31°7 M.G. 
30 B.D. + 23° 538 71 31 220% M.G. 
| 
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TaBLe II 


o’—a  cos(p—x) sin (p— x) cos (p—x) sin (p— x) 


” 


—2°2 +0°96 +0°28 + 1°00 —0°06 
— 0°25 +0°04 0°34 
+0°94 +0°34 + 0°63 —o'78 
+o°2 +0°98 +0°79 
+0°97 —0°23 + 0°99 +-O'17 
+0°99 + —o-98 
+o: +0°77 + 0°64 
—2'8 +o° —2°0 + 0°99 —O'l2 
+19 + 0° +0°27 + 0°96 
—2-7 + 0°99 —0'13 
+0'5 EF 1°00 + 0°02 
—o'8 +o" —o- + 1°00 
Lo» —2°2 +0'83 ~-0°56 


OVD DI DN 


The mean places for 1949-0 of the three stars (Nos. 2, 4 and 18) are given 
in Table III. For star Nos. 2 and 18, proper motions have been applied to the 
coordinates given in the catalogue in computing the mean place for use in the 
reduction. 


TABLE III 


R.A. (1949°0) Decl. (1949°0) Remarks 
h I , ” 


n ° 
23 46 7:97 — 443 91 A. G. Strasburg 8146 ; 
; ps0 (Yale) 
I 12 53°95 + 6 41 12°! Boss 1511 
4 21 36°79 +25 2 43°7 A. G. Berlin 1421 ; 
Hat 08-0027 ;45—0”-018 


Nizamiah Observatory, 
Hyderabad, Deccan, India : 
1950 February 4. 


ERRATUM 
M.N., 110, No. 1, 1950, P. A. Sweet’s paper: 


P. 76, equations (3.12) and (3.15), and the line immediately following equation (3.15), 
for w read W. 
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